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ABSTRACT 


The  MGD  generator  electrode  materials  problem  has 
been  greatly  simplified  by  Introducing  cooling  water 
into  the  electrode  stems.  Runs  of  over  one-half  hour 
duration  indicate  a  negligible  recession  of  the  surface 
so  that  the  same  electrodes  can  be  used  for  several  runs. 
Electrical  data  from  the  generator  Indicates  a  diffuse 
current  with  an  appj.led  electric  field  for  currents  up 
to  about  10  amperes  per  electrode  (current  density  of 
about  1-  amp/cm^).  Current  leakage  through  the  insula¬ 
ting  brick  has  been  examined  and  some  changes  ‘in  brick 
design  and  material  are  being  made  to  reduce  this  leakage. 

A  more  complete  report  than  that  previously  given 
is  made  of  the  investigation  of  tiie  Influence  of  the 
Ramsauer  effect  on  nonequilibrium  electron  temperatures. 
This  report  is  a  preprint  of  a  Research  Note  that  will 
appear  in  the  Physics  of  Fluids. 

The  calculation  of  electrical  conductivity  by  the 
Chapman -Enskog  method  has  been  compared  to  the  method  of 
Frost.  Very  satisfactory  agreement  has  been  shown  and 
the  assumptions  of  Frost  have  been  given  a  more  rigorous 
foundation. 

Diagnostic  methods  have  been  improved  with  the 
development  of  the  spectroscopic  line  reversal  technique. 
Temperature  measurements  believed  accurate  to  +  lO^K  have 
been  made  in  seeded  flames. 

An  investigation  of  the  effect  of  free  stream  tur¬ 
bulence  on  the  limiting  current  density  at  the  front 
stagnation  point  of  a  translating  cylinder  electrode  was 
undertaken.  An  effect  of  the  cylinder  diameter  and  whether 
the  electrode  was  cathodic  or  anodic  was  observed. 
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INTRODUCTION 


The  present  report  Is  the  fourth  Quarterly  technical 
report  written  under  supplemental  agreement  2(63-604)  to 
Contract  AF  49(638)~1123  which  is  supported  by  the  Advanced 
Reseai'ch  Projects  Agency.  Described  in  this  report  is  a 
research  effort  including  studies  of  magnetogasdynamlc  and 
fuel  cell  energy  conversion.  This  effort  involves  faculty 
and  students  in  the  departments  of  Mechanical  Engineering, 
Aeronautical  Engineering,  Chemical  Engineering  and  Chemistry. 
The  faculty  members  contilbutlng  to  the  program  are  as  follows. 


Professor  R.  H.  Eustls 
Professor  C.  H.  Kruger 
Dr.  M.  Mltchner 

Shock  Tube  Professor  D.  Bershader 

Fuel  Cells 

Kinetic  Studies  Professor  D.  Masor 

Dr.  P.  van  Rysselberghe 

Chemisorption  Studies  Professor  E.  Hutchinson 

In  addition  to  this  report,  tv;o  supplementary  reports 
are  included  which  conclude  a  particular  topic  of  research. 
These  reports  are  the  following: 

Supplementary  Report  Number  1,  "Heat  and  Mass  Transfer 
From  the  Surface  of  a  Cylinder  with  Discontinuous 
Boundary  Conditions  to  an  Incompressible  Laminar  Flow”, 
by  Zeev  Rotem  and  David  M.  Mason. 


Magnetogasdynamics 

Channel 


ill 


Supplementary  Report  Numler  2,  "Transport  Properties  of 

Partially  Ionized  Monatomlo  Oases",  Steve  deVoto  and 

Daniel  Bershader, 

with  this  report,  the  work  supported  by  oont’ act 
A?  49(638)-1123  on  fuel  cells  and  on  shock  heated  plasmas 
Is  concluded.  Research  In  these  areas  Is  summarized  in  the 
present  report  and  detailed  descriptions  are  given  In  the 
supplementary  reports  as  well  as  in  previous  Quarterly 
Reports.  Work  Is  continuing  on  studies  related  to  Magneto- 
sasdynamlos  energy  conversion  In  a  combustion  driven  generator 
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I. 


MAGNETOGASDYT'IAMICS 


A.  Continuous  Flow  Channel 


1)  Power  Generation  -  As  discussed  in  the  last 
progress  report,  the  most  satisfactory  overj.ll  generator 
pe.rformance  had  been  achieved  with  uncoated  type  316  stain¬ 
less  steel  and  molybdenum  electrodes.  The  estimated  surface 
temperatures  of  these  materials  during  operation  were  respec¬ 
tively  l800°F  and  2200°F,  and  the  exposed  surface  area  of  each 
electrode  was  0.4"  by  1-1/4". 

The  first  run  of  the  power  generation  apparatus  in  the 
current  report  period  incorporated  type  3l6  stainless  steel 
for  all  26  electrode  heads.  The  exposed  surface  area  of  each 
electrode  head  was  3/4"  in  the  flow  direction  by  1-1/4"  normal 
to  the  flow  direction.  Boron  nitride  Insulating  spacers  1/32" 
thick  were  used  between  the  backside  of  the  electrode  heads  and 
the  metallic  housing  to  increase  heat  tr-ansfer  and  reduce  the 
electrode  temperature.  In  addition  to  inves tigatl on  of  elec¬ 
trode  behavior,  the  purpose  of  this  test  was  to  measure  gas 
conductivity  and  leakage  resistances. 

The  following  nominal  run  conditions  were  used: 


combust ion -gas  flow  rate 
oxidizer  to  alcohol  ratio 
N^/Og  ratio 
seed  fraction 

inlet  stagnation  temperature 
static  pressure  at  exit 
gas  cross-sectional  area 
electrode  separation  distance 
magnetic  field  Intensity 


0.13  Ib/sec-in^ 
stoichiometric 
1.0  by  weight 
0.01  K  by  weight 
2700°K 
1  atm 

5-5  in^ 

4.4  in 
0 


Measurements  of  electrical  resistances  were  made  during 
the  preheat  and  pjst  heat  phases  (low  flow,  propane  fuel)  as 
well  as  during  the  principal  run.  Data  were  recorded  for 
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resistances  from  anode  to  ground  (test  section  housing ),  cathode 
to  ground,  and  anode  to  cathode,  in  all  13  circuits.  Electrical 
currents  during  the  preheat  phase  were  so  low  as  to  produce  no 
ammeter  deflections.  Indicating  leakage  resistances  of  at  least 
several  thousand  ohms  under  these  relatively  low  temperature, 
no-seed  conditions. 

The  first  12  minutes  of  the  principal  run  were  occupied 
by  the  usual  warmup  phase  fc  '*  lowed  by  recording  of  voltage- 
current  loading  characteristics  of  all  the  13  circuits.  The 
data  from  each  circuit  exhibited  linear  characteristics  up  to 
currents  of  8-10  amps  per  circuit.  The  single  data  point  for 
the  highest  applied  voltage  ■indicated  a  corresponding  lowering 
of  internal  resistance  at  10  of  the  13  circuits.  The  loading 
curve  for  the  central  circuit  is  shown  in  Figure  1.  The  average 
resultant  gas  resistance  at  each  circuit,  estimated  from  the 
slope  of  the  loading  curves  and  accounting  for  the  subsequently 
measured  leakage  reslotances,  is  approximately  24  ohm.  This 
resistance  corresponds  to  an  effective  gas  conductivity  of 
3.9  mho/meter.  Unfortunately,  the  conductivity  of  the  hot  gas 
core  cannot  be  determined  from  this  run  because  wiring  from  the 
sidewall-potential  probes  was  apparently  destroyed  early  in  the 
run  by  gas  flow  behind  the  sidev/all  brick. 

One  set  of  data  on  leakage  resistance  from  electrode  to 
ground  was  obtained  at  a  run  time  of  11. 5  min.  These  data, 
obtained  for  an  applied  potential  of  70  volts,  indicated 
corresponding  leakage  resistances  ranging  from  approximately 
30  oiim  to  400  ohm,  depending  on  circuit  location.  The  central 
circuits  had  the  highest  leakage  resistances.  Application  of 
a  higher  potential  between  electrodes  and  ground  resulted  in 
the  blowing  of  fuses  in  several  circuits  making  interpretation 
of  subsequent  data  difficult. 

Steady-state  temperature  data  from  t.hermocouples  In  the 
%0  sidewall  brick  were  extrapolat-d  to  give  a  surface  tem¬ 
perature  of  approximately  3500°?.  Eased  on  data  from  thermo¬ 
couples  inside  the  electrode  heads,  electrode  surface  temperature 
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ranged  from  approximately  1400®F  to  2000 depending  on 
electrode  location.  Total  nun  time  under  the  pnevlouely 
tabulated  flow  conditions  was  26  minutes. 

Examination  of  the  test  section  after  this  nun  showed  a 
dark  electrically-ccnductive  deposit  on  all  electrodes  and 
adjacent  MgO  crick  walls.  Recession  of  most  electrodes  was 
less  than  O.Ol”,  but  a  few  had  locally  receded  up  to  0.04". 
Several  pits,  apparently  caused  by  local  arcing,  were  present 
in  the  metallic  v/all  cf  the  test  section.  The  pits  seemed 
to  exist  at  gaps  between  the  boron  nitride  spacers  (see  above). 
In  all  cases  the  pits  were  adjacent  to  electrode  heads  In  those 
circuits  containing  t.he  fuses  blown  during  the  measurement  of 
leakage  resistances. 

Although  the  electrode  lifetime  in  this  test  was  adequate 
the  associated  deposition  of  contaralnant  on  the  inter-electrode 
brick  would  short  the  Hall  field  in  power  generation  runs. 
Thermal  analysis  shewed  that  stainless  steel  electrodes  cannot 
be  run  much  cooler  than  In  the  above  test  unless  the  cooling 
medium  is  brought  close  to  the  entire  heated  face  .  It  was 
decided  that  direct  water  cooling  of  the  electrodes  was  re¬ 
quired  . 

The  design  finally  selected  is  ill  strated  in  Figure  2. 

It  consists  of  a  molybdenum  head  cooled  by  condustion  to  a 
water-cooled  copper  sha.nk.  The  molybdenum  head  can  operate 
at  average  surface  temperatures  several  hundred  degrees  P 
lower  than  a  similar  stainless  steel  head,  due  to  the  large 
difference  in  thermal  conductivities  of  the  two  materials. 

The  head  and  shank  of  the  electrode  are  connected  by  a 
threaded  Joint;  powdered  graphite  is  added  to  this  Joint 
during  assembly  In  order  to  Improve  contact. 

The  subsequent  power  generation  run  incorporated  26  of 
the  electrodes  shovm  in  Figure  2,  The  head  size  was  again 
3/4"  by  1-1/4".  Two  types  of  molybdenum  were  used:  one  -was 
commercially  pure,  sintered  molybdenum  and  the  other  was  a 
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FIGURE  2  Water-Cooled  Electrode  Used  in  Run  of  8/19/64 


fcltanium-zlrconlum,  arc-cast  TZM  alloy.  MgO  spacers,  1/4" 
thick,  were  used  between  the  heads  and  the  test  section 
metallic  housing  (same  as  in  all  previous  tests  except  the 
last  one).  An  attempt  was  made  to  electrically  insulate  the 
test  section  and  diffuser  from  the  gas  flow  by  applying  a 
1000®?  strain  gage  cement  to  all  exposed  surfaces.  The 
diffuser  was  also  electrically  Isolated  from  the  test  section 
ground ) . 

The  following  nominal  conditions  were  used  in  this  run: 

combustion-gas  flow  rate 
oxidizer  to  alcohol  ratio 

^2/^2 

seed  fraction 

inlet  stagnation  temperaturr 
static  pressure  at  exit 
gas  cross-sectional  area 
electrode  separation  distance 
magnetic  field  intensity 

The  total  run  time  at  these  condltlcns  was  34  minutes 
including  a  6-1/2  minute  interval  during  vmich  the  seeded 
fuel  was  replaced  by  unseeded  fuel.  At  the  end  of  this  time, 
whe  run  was  autom.atlcally  stopped  by  a  flame-out  signal  from 
the  flame  sensor  which  was  caused  by  an  accumulation  of  opaque 
deposit  on  a  window  in  the  optical  path. 

The  electrodes  withstood  the  test  very  well,  better  than 
those  in  any  previous  test.  There  was  no  overall  recession 
which  had  been  typical  of  electrodes  in  previous  cests. 

Instead,  the  heated  faces  of  .most  electrode  heads  exhibited 
local  pitting;  this  pitting  occurred  over  the  entire  surface 
of  the  upstream  three  pair  of  electrodes.  There  was  frequent 
slight  head  damage,  giving  the  appearance  of  having  been  melted, 
in  localities  adjacent  to  the  hottest  regions  of  the  MjgO  brick. 
There  was  no  electrode  damage  at  the  interface  with  cooler 
regions  of  the  MgO.  The  heated  surface  of  a  few  electrodes 


0.16  Ib/sec-ln^ 
stoichlmetric 
1.0  by  weight 
0.01  K  by  weight 
2700®K 
1  atm 
t>.0  in^ 

4.0  in 

-  2 
2.0  weber/m 
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retained  a  foreign  material  which  appeared  to  have  flowed  from 
the  adjacent  upstream  MgO  spacer.  Also,  some  spacers  seemed 
to  have  experienced  a  buildup  of  material  during  the  test. 

There  was  no  deposit  on/or  melting  of  KgO  sidewall  brick. 
Comparison  of  molybdenum  -o  the  TZM  alloy  is  inconclusive. 

All  leakage  resistances  .measured  after  the  test  were 
several  craers  of  .magnitude  less  than  those  before  the  test. 

For  example,  the  indicated  resistance  between  two  adjacent 
electrodes  was:  before  test  -  2  Meg  Q  ;  after  test,  with 
inter-electrode  spacer  installed  -  1000  O  ;  after  test,  with 
inter-electrode  spacer  removed  -  4000  O  . 

Temperature  was  measured  at  one  internal  location  in 
each  of  three  molybdenum  heads,  two  in  upstream  circuit  #1 
and  one  in  downstream  circuit  #13.  Near  the  e.nd  of  the  run 
the  coi'respond ing  temperatures  were  approximately  1750*F, 

1350®F,  and  Temperatures  of  the  heated  surfaces  were 

probably  within  100®F  of  these  recorded  temperatures. 

Internal  temperatures  were  also  measured  of  upstream 
MgO  brick  in  the  test  section.  These  te.mperatures  were 
relatively  constant  after  warmup.  Extrapolations  indicate 
surface  temperatures  of  up  to  2800‘’F  on  the  inter-electrode 
spacers  and  about  3500®F  on  the  sidewall. 

Figures  3  through  10  present  electrical  data  showing 
representative  behavioi  of  the  13  separate  lead  circuits. 

More  detailed  information  is  presented  for  two  circuits, 

#2  and  #12;  these  circuits  typify  the  general  behavior  of 
circuits  in  the  two  halves  of  the  generator.  The  open-circuit 
load-potential  histories  shown  in  Figure  3  are  consistent  with 
previous  results.  Data  indicate  that  as  time  Increases  a  gradual 
decrease  in  voltage  occurs  for  circuits  in  the  upstream  1/2  to 
2/3  of  the  generator  but  the  voltages  remain  relatively  constant 
in  the  downstream  circuits.  The  short-circuit  current  histories 
shown  in  Figure  4  are  not  as  steady  as  exu'^ctei.  As  in  pr^-vious 
tests,  approximately  six  minutes  were  requires  to  reach  t.ne 
apparent  steady  level;  this  is  consistent  with  temperature 
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data  from  the  sidewall  brick  and  cooling  water.  Comparison  of 
electrical  data  obtained  before  and  after  operation  with  un¬ 
seeded  fuel  have  lead  only  to  the  conclusion  that  if  a  change 
in  leakage  resistance  occurred  during  the  6-1/2  minute  Interval 
then  its  effect  on  open  circuit  voltage  was  limited  to  approxi¬ 
mately  2  volts. 

Figure  5  presents  the  axial  distribution  of  open-circuit 
voltage  available  to  a  load.  At  early  times  shorting  apparently 
occurred  at  the  center  of  the  generator.  All  voltage  distribu¬ 
tions  after  3  minutes  of  running  were  similar  to  that  shown  for 
31.5  minutes.  Axial  distributions  of  load  current  at  the  four 
loads  used  In  this  run  are  shown  in  Figure  6.  An  axial  distri¬ 
bution  of  electrode  potential  relative  to  the  grounded  test 
section  housing  Is  shown  in  Figure  7  for  operation  with  all 
circuts  essentially  short  circuited.  Figures  5  and  and 
perhaps  Flgsre  6  demonstrate  the  presence  of  local  shcrsing 
within  the  generator.  The  difference  between  electrode  poten¬ 
tials  at  circuts  #13  and  #1  is  herein  called  the  total  Hall 
potential.  The  measured  history  of  the  total  Hall  potential 
for  those  runs  with  all  load  circuits  shorted  is  shown  in  Figure 
8.  These  were  the  most  unsteady  data  observed  In  the  entire  3^ 
minute  run.  The  expected  total  Hall  potential  is  300  volts, 
based  total  Hall  potential  is  300  volts,  based  .n  predicted 
gas  properties  and  measured  currents.  From  Figure  7  it  is  seen 
that  the  expected  Kail  field  was  actually  achieved  between  elec¬ 
trodes  #7  and  #8  and  also  befween  electrodes  #11  and  #13. 

Figures  9  and  10  show  the  loading  characteristics  of  the 
two  circuits  for  which  histories  were  shown.  Most  circuits 
produced  more  repeatable  data  and  some  .had  significantly  more 
linear  characteristics  than  those  shown. 

Comparison  of  loading  characteristics  and  current  distri¬ 
bution  for  all  13  circuits  has  indicated  a  sort  of  equilibrating 
effect  among  adjacent  electrodes.  Frequencly  a  circuit  which 
performed  less  than  the  tre.nd  was  adjacent  to  a  circuit:  which 
performed  better  than  the  trend.  It  is  therefore  meaningful 
to  investigate  the  cumulative  behavior  of  ail  circuits  in  the 
generator.  This  is  done  in  Figures  11  and  12. 
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Circuit 

C3 
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FIGURE  10  Loading  Characteristics, 


Circuit  #12 


Figure  11  presents  the  averaged  loading  curve  for  the 
total  generator.  The  "Averaged  Load  Potential"  is  simply  a 
numerical  average  of  all  the  13  simultaneous  load  voltages, 
and  the  "Total  Load  Current"  is  the  sum  of  the  13  simultaneous 
load  currents.  The  straight-line  correlation  of  the  data  is 
indicative  of  constant  generator  internal  resistance,  inde¬ 
pendent  of  lOv^i  current .  Figure  12  shows  the  corresponding 
total  generator  power  output  at  the  various  loads.  The 
prediction  is  for  a  simple,  constant  internal  resistance 
generator,  namely 


Power  =  V 


oc 


Ave 


■•SC 


SC 


The  average  open  circuit  load  potential,  and  the 

total  short  circuit  current,  I  ,  are 

165  volts  and  42.8  amps,  as  obtained  from  Figure  11.  It  is 
noted  that  the  above  form  of  the  equation  for  generator  power 
is  independent  of  the  degree  of  electrode  segmentation. 

Figure  13  presents  typical  gas  pressure  data  at  open 
circuit  and  short  circuit. 

A  few  measurements  of  end  effects  were  made  by  short- 
circuiting  all  Interior  electrodes  ana  successively  open- 
circuiting  the  exterior  electrodes.  The  results  showed  a 
significant  increase  in  load  current  for  the  two  end  electrodes 
of  the  Interior  set.  The  number  of  interior  electrodes  which 
were  affected  increased  as  the  operating  length  of  the  generator 
was  reduced.  For  exam.ple,  with  circuits  #1  through  #4  and  #10 
tnrough  #13  open-circuited,  t.he  active  generator  consisted  of 
circuits  #5  through  #9  operating  in  the  short-circuit  condition. 
For  this  case  the  current  through  circuit  #5  (upstream  end) 
was  30^  higher  than  when  a^  13  circuits  were  shorted;  the 
corresponding  increase  for  circuit  #7  (center)  was  11%  and 
for  circuit  #9  (downstream  end)  was  16,^. 

Several  measurements  of  generator  internal  leakage 
resistances  were  made  during  the  6-1^ ^2  minute  Interval  of 
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Short  Clr 


FIGURE  13  Combustion-Gas  Pressure  X-Dlstributlon 


operation  with  unseeded  fuel.  Qualitative  measurements  were 
made  of  resistance  from  electrode  to  ground,  between  adjacent 
electrodes,  and  between  opposing  electrodes.  Indicated 
resistance  varied  from  30  ohm  to  600  ohm.  All  measurements 
were  made  In  the  vicinity  of  the  electrodes  In  circuit  #7. 
Although  the  measurements  should  not  be  used  quantitatively, 
they  Indicate  leakage  resistances  of  the  same  order  of 
magnitude  as  the  overall  generator  resistance  and  the  gas 
resistance.  Slopes  of  the  loading  characteristics  curves 
such  as  in  Figures  9^  10,  and  11  Indicate  an  overall  generator 
average  resistance  of  approximately  50  ohm  per  circuit. 

The  potential  between  test  section  and  diffuser  was 
measured  at  several  times  during  the  run.  In  ail  cases,  no 
potential  was  measured.  Much  of  the  ceramic  coating  was 
removed  from  the  diffuser  wall  during  the  run.  This  apparently 
had  relatively  little  effect  on  the  generator  performance 
because  the  expected  Hall  field  seemed  to  be  developed  In  the 
downstream  electrode  region. 

Very  limited  data  were  obtained  with  the  electrodes, 
connected  for  continuous  electrode  operation.  Automatic 
shutdown  occurred  soon  after  the  electrical  ccn.nectlons  were 
made.  However,  the  open-circuit  data  clearly  show  the  effect 
of  the  Internal  Hall  current.  .The  data  also  indicate  a  low 
resistance  path  to  ground  in  the  vicinity  of  electrode  #10; 
this  is  consistent  with  low  load-current  in  circuit  #10  (see 
Figure  6). 

The  next  run  in  the  power  generation  program  will 
concentrate  on  the  leakage  resistance  problem.  High  density 
low  porosity  Al20^  will  be  used  in  several  locations  In  order 
to  reduce  absorption  of  conductive  materials  by  the  brick. 
Various  Joints  will  be  tried.  The  electrodes  from  the  last 
run  will  be  re-used  with  the  threaded  Joint  between  head  and 
shank  replaced  by  a  brazed  Joint. 
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2)  Heat  Transfer  and  Fluid  Friction  -  Three  problems 
are  being  considered  in  the  area  of  heat  transfer.  The  first 
is  cone erne,  with  the  insulator  boundary  layer  in  an  MGD 
channel  and  the  second  treats  the  electrode  boundary  layer. 
These  analyses  take  into  account  details  of  the  MGD  inter¬ 
action  and  in  their  early  stages  ignore  affects  such  as 
compressibility  or  chemical  reaction.  A  third  study  is 
being  made  of  the  heat  reansfer  to  a  cooled  electrode  or  heat 
transfer  measuring  calorimeter  and  is  principally  concerned 
with  compressibility  and  effects  due  to  large  temperature 
differences  between  the  gas  and  wall. 

The  proceeding  progress  report  gave  a  detailed  discussion 
of  this  ''esearch  and  this  work  has  continued  actively  during 
the  present  quarter.  The  new  results  from  these  investigations 
however,  are  incomplete  and  a  detailed  account  will  be  given  in 
the  following  Quarterly  Report. 

3)  Physics 

a)  Calculation  of  Transport  pi^ODerties  in  Seeded 
Combustion  Gases  -  In  I96I,  Frost^,  propos-d  calculating  the 
electrical  conductivity  of  a  partially  ionized  gas  using  Allis’ 
equation,  valid  for  a  Lorentzian  gas.  Frost  replaced  the 
electron-neutral  collision  frequency  in  Allis'  formula  with 
the  sura  of  the  electron-neutral  and  electron-ion  collision 
frequencies.  To  take  into  account  electron-electron  inter¬ 
actions,  Frost  modified  the  expression  for  the  electron-ion 
collision  frequency  sc  as  to  obtain  the  Spitzer-Harm  value 
of  the  conductivity  in  the  limiting  case  of  a  fully  ionized 
gas. 

In  order  to  examine  the  validity  of  Frost’s  procedure, 
we  have  determined  the  electrical  conductivity  of  a  partially 
ionized  gas  taking  into  account  rigorously  e-e  and  e-l 
collisions  by  a  modification  of  the  Chapman-Enskog  expansion 
scheme . 
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For  the  simple  model  of  electron-neutral  collision 
frequency 

m  m 

''en  =  'V  (i-)  '  (1) 

and  the  classical  electron-alectrcn  and  electron-lcn  cross- 
sections,  the  first  three  approximations  to  the  electrical 
conductivity  have  been  derived  for  an  arbitrary  degree  of 
Ionization.  All  three  approximations  agree,  with  the  knov.n 
Corresponding  approximations  in  the  two  limits  of  either  a 
Lorentzian  or  a  fully  ionized  gas. 

Tn  determinant  form  the  third  approximation  may  be 
v.Tltten  as: 
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The  quantities  and  o_,  are  the  exact  conductivities 
of  the  lA-‘rentzian  gas  and  the  fully  ionized  gas  respec^i.el 
The  remaining  quantities  depend  on  the  assumed  power  lav;  of 
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Figures I4  to  I7  show  the  dependence  on  degree  of  ionl 
of  the  first  three  approximations  to  the  electrical  cond. 
These  curves  give  a  clear  idea  of  the  rate  of  convergence 
first  three  approximations  for  different  models  of  the  el 
neutral  cross-section.  For  high  values  of  ionization,  el 
interactions  dominate.  The  conductivity  is  insensitive  t 
parameter  m  and  rapidly  approaches  the  value  of  the  cor 
ponding  fully  ionized  case.  The  behavior  at  low  levels  0 
ionization  is  different  however.  As  shown  by  the  detaile 
study  of  the  Lorentzlan  gas  in  the  last  progress  report^, 
the  speed  of  convergence  is  .hignly  sensitive  to  the  depen: 
of  cross-section  on  lectron  velocity.  For  m  =  -  1  ar.''* 

m  =  0  (Maxwellian  molecules)  convergence  is  excellent  and 
third  approximation  seems  to  be  sufficient  for  all  levels 
ionization.  For  m  =  1/2  (hard  spheres)  and  for  m  =  1 
convergence  for  low  values  of  ionization  is  less  rapid. 
While  the  third  approximation,  for  these  two  rodels,  is 
certainly  satisfactory  at  high  levels  of  ionizations,  it  i 
not  very  good  at  low  ionization  levels. 

At  low  values  of  ionization  Frost's  results  seem  to  c 
more  satisfactory,  as  expected,  since  his  formulation  is 
rigorous  for  low  ionization  and  corresponds,  in  the 
to  the  infinite"  Chapman-Enskog  approximation,  o(«). 
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FIGURE  14  Normalized  Conductivity  of  Pii-at  Three  Chapman -En 3 kog 

Approximations  and  of  Frost  as  a  Function  of  Ionization. 
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FIGURE  lo  Normalized  Conductivity  of  First  Three  Chapman -Ena kog 

Appro;: imat Iona  and  of  Frost  as  a  Function  of  Ionization. 


Normalized  Conductivity  of  First  Three  Chapman -Ensk'og 
Approximations  and  of  Frost  as  a  Function  of  Ionization 


For  high  values  of  ionization.  Frost’s  results  agree  well 
with  the  third  Chapnian-Enskog  approximation  which  includes  the 
electron-electron  interaction  rigorously.  Since  Frost  has 
used  a  fitting  tf-hnique  to  make  his  results  agree  with  the 
values  of  the  fully  ionized  conductivity,  his  results  for 
infinite  degree  of  ionization  are,  as  expected,  slightly 
(less  than  2%)  better  than  t.he  third  approximation. 

In  conclusion,  we  believe  that  this  study  provides 
thoretical  Justification  for  the  use  of  Frost's  method  in 
calculating  the  electrical  conductivity  of  a  partially 
ionized  gas  in  the  absence  of  a  magnetic  field. 

A  second  related  problem  was  investigated  -  that  of  the 
effect  of  small  in  a  .  T.he  kinetic  theory  description  of 
collisions  between  charged  particles  has,  until  recently^, 
required  some  sort  of  ad  hoc  cut-off  procedure  in  order  to 
prevent  t.he  occurrence  of  divergent  integrals.  The  source 
of  this  difficulty  stems  from  the  long-range  nature  of  the 
Coulomb  Interactlo.n  ana  the  treatment  of  all  collisions  as 
two-body  encounters.  In  actuality,  the  interaction  potential 
is  essentially  shielded  by  t.he  collective  behavior  of  the 
particles  at  sufficiertly  large  distances.  A  convergent 
theory  is  obtained  by  i.ntroducing  some  device  for  ignoring 
collisions  wi^h  large  i.mpact  para.meters  in  excess  of  the 
Debye  lengt.h.  T.he  resulti.ng  theory  is  valid  to  order 
(in  h)”^  ,  where  h  =  1.24  x  10^ is  the  ratio  of 
the  Debye  length  to  t.he  impact  parameter  for  a  90°  deflection. 
(The  quantities  T  and  n  aenute  the  temperature  in  °K  and 
the  number  density  per  cm^  respectively).  In  particular,  the 
widely-quoted  value  of  Spitzer  and  Har.m  for  the  electrical 
conductivity  of  a  fully-ionized  pla.sma  is  correct  t 
order. 


t.hl: 


In  the  folluwing  table,  we  have  calculated  the  values  of 
A  and  in  a  for  a  range  of  te.mperatures  a.nd  number  derisities 
whlcn  en  ^npass  conditions  expected  in  foHD  generators.  Typical 
values  of  in  .v  are  seen  to  range  bet’ween  s  and  5  which  Imolie 


are  seen  to  range  between  4  arid  5 


an  uncertainty  of  about  25^  In  the  Spltzer-Harm 


conductivity 


Values  of  A  and  ( /n  a> 


ng(cm"^) 


t^Ck) 

1  10^'^ 

10^^ 

loi'* 

1015 

124 

(4.82) 

39.2 

(3.67) 

12.4 

(2.52) 

3.92 

( 1 .  j7  ) 

— 

72. 1 

(4.28) 

22.4 

(3.11) 

7.21 

(1.98) 

2000 

1110 

(7.01) 

111 

(4.71) 

35.1 

( J06) 

1= - 

11.1 

(2.40) 

2b00 

1550 

(7.35) 

491 

(6.20) 

155 

(5.04) 

49 . 1 
( 3 . 69 ) 

15.5 

(2.74) 

3000 

2040 

(7.62) 

645 

(6.47) 

204 

(5.32) 

64.5 

(4.17) 

20.4 

(3.02) 

The  region  to  the  left  of  the  heavy  line  In  this  table 
ndloatea  roughly  conditions  for  ..hlch  this  theory  la  auffl- 
clently  accurate  (uncertainty  less  than  20^). 

as  „elf  Rosshbluth 

ically  took  into  account  the  collective  behavior  of  a  plas-na 

and  converged  for  large  Inpact  parameters.  However  thL 

eo-called  "fluld-appro.lmntlon"  dlv -rged  for  amaU  C 

par^eters,  and  so  It  was  necessary  to  Introduce  an  ad  ti  o 
imll  cut-off  limit. 

of  ch  Important  progress  In  the  klnetlc-theorj 

Charged  particles  was  made  In  196j  by  Klhara  and  Aono'* 

LTllL"  :  ^  -"^i-ng  b^th 

e  and  distant  encounter  contributions  to  the  kinetic 
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descz’iption  of  collisions,  in  such  a  way  as  to  yield  a  diver¬ 
genceless  theory  without  the  need  of  ad-hoc  assumptions.  The 
validity  of  this  new  theory  is  of  order  ,  and  thus 

provides  ample  accuracy  for  MHD  generator  applications.  The 
region  to  the  left  of  the  double-line  in  the  preceding  table 
indicates  conditions  for  which  this  new  theory  Is  sufficiently 
accurate  (uncertainty  10^  or  less). 


The  theory  of  Klhara  and  Aono  has  been  applied  by 
Ikitawa^  to  the  calculation  of  the  electrical  conductivity 
of  a  fully-ionized  plasma.  The  method  of  calculation  is 
similar  to  that  of  Chapman-Enskog,  and  involves  an  expansion 
in  terms  of  Sonlne  polynomials.  In  Figure  i8  we  have  plotted 
the  ratio  of  Ikltawa's  second  approximation  to  the  Spitzer- 
Harm  value  of  the  electrical  conductivity,  as  a  function  of 
in  A  .  For  vaiues  of  in  A  4  ,  the  more  accurate  value  of 
the  electrical  conductivity  is  about  3C^  higher  than  the 
Spitzer-Harm  value . 


In  1964,  Hynn^  reported  a  measurement  of  the  electrical 
conductivity  in  a  fully-ionized  plasma  whicn  was  20^  higher 
than  the  Sp-Ltsspr-Harm  theoretical  value.  The  probable  error 
or  Hynn‘3  measurement  was  stated  to  be  +  10«.  At  the  value 
of  in  =  7  fsr  which  Rynn’s  experiment  was  performed, 
Ikitawa's  value  of  conductivity  is  about  lOp  larger  than  the 
Spltzer-Karm  value,  which  places  Ikitavfa‘s  theoretical  value 
within  the  prcbabls  error  of  the  experiment.  This  comparison 
with  experiment  appears  to  provide  some  support  for  . the  new 


theory. 

Under  MiD  generator  conditions,  the  plasma  is  only 
partially  ionized,  and  so  the  magnitude  of  the  correction 
discussed  above  will  be  proportionately  smaller.  When  the 
charged  particle  encounters  are  equally  Important  to  electron- 
neutral  collisions,  we  may  expect  a  correction  of  the  order 
of  15^,  and  so  we  plan  to  incorporate  the  results  of  this  new 
theory  into  our  calculations. 
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b)  The  Influence  of  the  Ramsauer 

Electron  Temperatures  -  The  following  is  a  Research 
Note  submitted  for  publication  in  the  Physics  of  Fluids.  We 
note  that  recent  calculations  by  Zukoski,  Cool,  and  Gibson 
(AIAA  Journal  1410)  confirm  the  important  influence  of  the 


Ramsauer  cross-section, 


,7,8/9 


Kerrebrock.  and  others*  have  considered  recently  the 
effect  of  elevated  nonequllibrium  electron  temperatures  on  the 
electrical  conductivity  of  partially-ionized  seeded  noble-gas 
plasmas.  The  departure  of  the  electron  temperature  from  the 
gas  temperature  has  been  estimated  by  means  of  a  balance 
between  the  energy  gained  by  the  electrons  as  they  drift 
through  the  electric  field  and  the  energy  lost  in  collisions 
with  the  heavy  particles.  Constant  collision  cross-sections 
have  been  used  in  these  calculations  and  in  the  calculation  of 
the  electrical  conducoivity .  The  electron  number  density  has 
been  found  by  mea.ns  of  the  Saha  equation  evaluated  at  the 
electron  temperature  and  this  has  been  Justified  by  the 
assumption  of  a  Maxv/ellian  dis tr:.  jution  of  electron  energies. 
In  the  present  analysis  a  displaced-Maxv/ellian  electron  dis¬ 
tribution  functlcn  is  used  to  assess  the  imoortance  of  the 
variation  of  the  cbilision  cross-section  with  electron  energy. 
In  particular,  for  elastic  collisions  between  electrons  and 
noble -gas  atoms  the  pronounced  minimum  in  the  cross-section 
at  electron  energies  of  the  order  of  one  electron  volt  (the 
Ramsauer  effect)  is  shown  to  have  a  significant  effect  upon 
the  variation  of  the  collision  frequency  with  electron  tem¬ 
perature  . 

We  adopt  a  coordinate  system  for  which  the  mean  velocity 
of  the  heavy  particles  is  zero  and  assume  a  Maxwellian  heavy- 
particle  distribution  function  at  the  known  temperature  T^  . 
The  electron  distribution  function  is  taken  to  be  Maxellian 


about  the  electron  drift  velocity  c  with  an  electron 

e 


temperature  T 


The  unknown  parameters  c^  and  T^  are 


found  by  means  of  the  momentum  and  energy  moments  of  the 


-electron  Boltamann  equation  with  only  the  eleotrlo-fleld  term 
=3P#tained  on  the  left-hand  side. 

With  the  displaced-Maxwellian  distribution  our  collision 
Integrals  are  identical  to  those  considered  by  Morse^^.  For 
small  values  of  the  parameter  ni^cp(2k  T„)  and  with  the 
definitions  J  s  -  em^c^  and  a  s  J/B  ,  the  results  can  be 

put  in  a  form  that  is  similar  to  familiar  phenomenological 
relations.  We  then  have 


“e  4 
e^7  3 


2  n^  Q  ,  (T  ) 
h  ^e,h  ^*e' 


o  e 


3k(T 


■  ^h)  I  J 


where  the  effective  mean  cross-section  (t  )  is  defined 
In  terms  of  the  differential  cross-section^  n^S^X)  by 


;  (T  .3  -  -m^g  /2kTg  ^ 

§ih  e  2icTg  ^  ^  S  2w  /(l-cosx}lg^j^(g,X}sinxdxdg  . 


%,h^*^e^  differs  by  only  a  factor  from  the  integral  (l) 

that  enters  tne  Chapman -Enskog  theory.  The  summation  is^r- 
^ed.  over  all  species  of  heavy  particles  (neutrals  and  ions), 
W  gLr^^^^  particles  are  equal,  gqs.  (3)  and 


^  3K(T  -  TJ 
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The  effective  cross-section  evaluted  for 

several  noble  gases ^  potassium,  and  cesium  by  numerical  inte¬ 
gration  over  the  relative  velocity  g  .  Collision  probability 
data  were  taken  from  Brown"^  and  Massey  and  Burhop  and  the 
references  cited  therein.  Prom  estimates  based  on  the  avail¬ 
able  data  for  the  angular  dependence  of  the  differential 


cross-section,  the  momentum- transfer  cross-section 

2Tr  (l*'COSx)l^  (g,x)sinxdx  in  Eq.  (5)  «as  approximated 

by  the  measured  total  cross-section  2tr  /  I  * 

Our  calculations  indicate  ttiat  the  effect  of  uncertainties 

in  the  very- low-energy  cross-section  is  snail  for  T  >2000®K; 

this  results  from  the  factor  g"  in  Eq.  (p)*  Tne  calculated 

values  of  Q  given  in  Figure  19-  gases  that 

exhibit  a  pronounced  Ramsauer  effect  the  effective  cross- 

section  Q  ^  varies  by  as  much  as  a  factor  of  8(argon)  in 

the  range  of  electron  temperatures  shwwn.  These  values  of 

Q  .  are  significantly  greater  than  the  oocal  cross-section 
0  ^ 

at  the  Ramsauer  minimum.  The  effective  eiectron-ion  cross- 
section,  calculated  from  Eq.  (5)  with  the  Rutherford  cross- 
section  and  the  Debye  cutoff,  is 


(Tg) 


2(kT^)^ 


(7) 


With  this  displaced-Maxwelllan  cross-section  the  electron-ion 
resistivity  (Eq.  3)  is  higher  than  the  Soitzer-Harm  resisGivloy 
by  a  factor  of  1.98. 

Calculations  based  on  Eqs.  3j  5  s.n6  6  and  che  Saha 
equation  evaluated  at  the  electron  temperature  were  compared 
with  the  measurements  of  Kerrebrock^  and  Zukoskl'^.  These 
experiments  were  conducted  with  flowing,  partially-ionized, 
potassium-seeded  argon  at  atmospheric  pi'essure  with  an  applied 
electric  field.  The  comparisons,  shown  in  Figures  20  and  21, 
are  favorable  in  view  of  the  fact  that  no  corrections  to  the 
theory  or  empirical  information,  other  than  measured  cross 
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EISCTRON  temperature,  Tg  ,  In  10^  "K 

FIOURE  19  Calculated  effective  cross-aectlons 
a  i’unotion  ol  electron  temperature. 


ELECTRICAL  CONDUCTIVITY,  a  ,  mho/cm 


10 


CURRENT  DENSITY,  J  ,  amps/cm^ 


FIGURE  20  Comparison  of  the  present  theory  ( - ) 

with  the  experiments  (■ - -•)  of  Kerrebrock^. 

T^  ss  2000  *’K. 
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FIGURE  21  Comparlsou  of  the  present  theory  (-- . . )  with 

the  experiments  ( - )  of  Zukoakl^.  Tk  *»  2000®K 


sections,  were  used  .  Calculation,  made  in  reference  8 


nd 


9  with  the  constant  cross -sections 


■16 


and 


Q  -  0.7  X  ^0 
-lb  2  " 

=  250  X  10  cm  lie  slgnif icantly  above  the  data. 
Although  our  results  indicate  the  Importance  of  the  variation 
of  the  differential  cross-section  with  electron  energy  on 
nonequilibrium  electron  temperatures  and  tend  to  support 
use  of  a  displaced-Maxwelllan  distribution,  the  situatic 
not  be  fully  understood  until  the  electron  distribution 
tion  can  be  accurately  predicted.  Such  an  investigatior; 
currently  in  progress. 
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4)  Diagnostics  -  A  sodium  line  reversal  apparatus 
has  been  developed  in  conjunction  with  the  Stanford  MGD 
generator  channel  to  provide  a  check  on  the  calculated 
temperature  and  to  measure  directly  the  electron  temperature 
in  the  presence  of  induced  currents.  A  description  of  the 
method  and  its  application  to  the  Stanford  channel  is  given 
in  detail  in  Appendix  A  of  this  report. 

A  llg.ht  source  which  is  a  calibrated  tungsten  ribbon  lamp 
is  used  in  conjunction  with  an  optical  system  and  a  Jarre  11- 
Ash  half  meter  Ebert  scanning  spectrometer.  Experimental 
results  obtained  wit.n  the  present  apparatus  have  been  per¬ 
formed  on  small  prepane-air  torch  arrangements,  with  small 


amounts  of  sodium  added 


by  means  of  a 


rod  coated  with  NaD  inserted  into  the  base  of  the  flame. 
In  this  phase  of  development  the  temperature  of  the  flax 
has  been  determined  to  an  observable  experimental  accur.; 
of  10®K.  By. this  it  is  meant  that  two  different  output 
spectral  line  shapes  corresponding  to  a  difference  of  1C 
in  flame  temperature  are  easily  distinguished  and  accura 
measured.  Figures  A-2,  A-3,  A-4,  and  A-5  in  Appendix  A- 


'  y-!.  ^ 

yO  Oi 


The  fact  that  losses  by  radiation  and  the  effec' 
collisions  were  not  included  in  the  analysis  must  be  cc 
however,  in  the  interpretation  of  these  comparisons. 


—  ^ 

°K 

tely 
1  are 

^  1  p  ^  4 
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topical  examples  of  the  NaD  lines  as  seen  in  the  small  propane 
flame , 

Attempts  v/ere  made  to  obtain  line  reversal  with  a  potas¬ 
sium  seeded  small  propane  flame  for  the  potassium  lines  at 

O 

4042  and  4044  A  .  The  signal  received  from  these  lines  was 
very  weak  due  to  the  small  amounts  of  potassium  added  which 
the  physical  size  of  the  flame  dictated.  Emission  lines  were 
obtained,  but  the  reversal  condition  signal  was  immeasurable. 

At  the  present  time  testing  is  underway  to  determine 
flame  temperature  for  various  constituent  percentages  in  the 
MHD  generator.  The  temperatures  are  obtained  at  the  plenum 
chamber  exit,  exhausting  to  the  atmosphere.  The  preliminary 
results  from  this  testing  show  that  the  percentage  of  sodium 
in  the  flame  should  be  in  the  order  of  1/100  percent  or  less 
to  avoid  self  absorption  in  the  emission  lines.  Reversal 
conditions  have  been  observed  for  several  constituent  per¬ 
centages  at  low  flow  rates. 


B.  SLock  Tube 


■  The  folloving  paragraphs  constitute  a  brief  sunraary  of  the 
shock  tube  phase  of  the  ^5HD  energy  conversion  p-rogram. 

The  function  of  the  shock  tube  work  has  been  to  perform  quanti¬ 
tative  experiments  with  shock-heated  plasmas  which  would  axigment  our 
understanding  of  their  physical  properties  and  dynamic  behavior.  The 
tube  used  her^  was  designed  to  cover  and  extend  somewhat  the  range  of 
parameters  having  importance  for  the  performance  of  MHD  power  generation 
equipi!«nt.  Details  of  the  design  have  been  given  in  previous  reports  of 
this  series.  Performance  of  the  device  h.as  been  in  close  accord  with  its 
design,  and  has,  in  general,  been  very  satisfactory.  It  remains  the  most 
quantitative  device  available  for  producing  a  well-defined  sample  of 
plasma  having  precise  physical  and  dynamic  characteristics  in  the  tem¬ 
perature  and  density  range  of  interest  to  this  project. 

The  major  initial  effort  in  the  program  concerned  the  above- 
mentioned  shock-tube  design,  acquisition  and  fabrication  of  components, 
sheikedown  tests  of  subsystems,  and  a  series  of  calibration  shots. 

Initial  operation  took  place  in  mid-fall  1962  according  to  sched’ole,  and 
the  facility  has  seen  regular  use  since  that  time  in  a  series  of  research 
investigations.  By  now,  the  following  results  have  been  accomplished: 

.  (1)  Calibration  of  shock  tube  performance  over  a  wide  range 

of  operating  conditions; 

(2)  Evaluation  of  the  overall  precision  and  limitations  on 

the  use  of  optical  interferometry  to  measure  electron 

15  l8 

density  in  the  dense  plasma  range,  10  ^  -  10  /cc  ; 

(3)  Measurement  of  the  thennal  equilibration  of  a  shock- 
heated  plasma; 

(k)  Calculation  of  the  high  temperature  transport  properties 
of  moriatomic  gases,  using  the  higher  approximations  of 
the  Chapman-Enskog  forroaJLism. 

In  addition,  a  series  of  systematic  investigations  is  now 
under  way  in  the  following  areas: 
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1  ^ccaflJiisatlon  processes  in  a  floving  plasisa  cooled  by  a 
;.^8tatlQnafY  expansion  wave; 

-  -  C^)  fediative  thersal  timnsport. 

-  portions  of  the  above  stnnloe  have  been  given  In 

«rlltr  reports  in  this  series.  In  particular,  itens  (2)  and  (3) 
were  treated  in  Supplement  2  to  the  Qtmrterly  Report  1  March  -  31  May 

below).  In  particular,  this  work  did  verify  the  high 
preeiaion  of  optical  Interferometry  for  a»asuring  the  electron  density 
of  dense  plasmas.  It  also  verified  the  atomic  collision  model  leading 
to  theiaal  ionisation,  with  indication  of  the  order  of  mgnitude  of  atom- 
at^  «>llision  cross  sections.  The  work  also  yielded  results  for  the 

^«f^rature  T^rlation  in  tne  non-equilibrium  regime,  and  showed 
at  what  stage  of  the  equilibration  process  recombination  events  become 
l^rtant.  The  variation  of  the  electron  density  downstream  of  the  equi- 
Ut^tion  point  gives  further  quantitative  data  on  the  radiation  cooling 
rate  in  an  otherwise  uniform  plasna. 

described  also  provided  Implicit  verification  of 
-  equation,  in  that  the  measured  equilibrium  value  of  electron 

density  was  in  good  agreement  with  that  predicted  by  a  combination  of 

^nd  Saha  relations  applied  to  the  measured  shock 
and  shock  tube  conditions.  This  result  was  fiirther  verified  by 
an  additiot^  liBiependent  series  of  experiments  which  formed  the  basis 
of  aa>ther_ student’ s  Rngineer’s  ^nesis  (see  listing  below). 

The  analysis  of  high  teugjerature  transport  prepertles  of  a 
partially  ionized  monatomic  gas,  item  (4)  above,  has  helped  to  resolve 
the  disagreea^nt  which  develops  when  partially-ionized  conductivity 
ej^Jresslotts  are  extended,  in  the  limit,  to  full  ionization  -  the  trouble 
being  that  these  do  not  agree  with  results  calculated  explicitly  for  a 
fully  ionized  plasm.  The  present  work  has  shown  that  satisfactory 
agreement  can  be  obtained  If  one  goes  to  the  third  approximtion  in  the 
Chapaan-lnskog  tMory.  depressions  have  been  obtained  to  the  fourth 
^roxlaation  in  the  thex^  conductivity,  thermal  diffusion  coefficient 
aM  ordinary  diffusion  coefficient.  The  viscosity  Is  treated  to  the 

a^rox^tlOT.  It  Is  noteworthy  that  thermal  diffusion  plays  an 


important  role  in  an  ionized  gas. 

In  connection  with  Item  (5)  above,  recent  interferograms  have 
shown  the  need  for  obtaining  greater  fringe  shifts,  since  it  is  desired 
to  measure  the  variation  of  fringe  shift  along  a  radial  line  of  the 
Prandtl -Meyer  expansion  and  along  the  streamline  running  close  and  paral¬ 
lel  to  the  downstream  wall.  Ti'.e  total  fringe  shift  plot  given  in  the 

ft 

last  quarterly  progress  report  was  for  a  wavelength  of  45C0A  .  This 

O  O 

has  now  been  extended  to  wavelengths  of  5390A  and  69OOA  (correspon¬ 
ding  to  the  radiation  from  a  ruby  laser).  Since  the  major  contribution 
to  the  fringe  shift  is  from  the  free  electrons,  it  turns  out  that  tne 
higher  wavelengths  yield  greater  fringe  shifts  (proportional  to  weve- 
length) . 

Turning  to  item  (6),  we  remark  only  briefly  that  preliminary 
experiments  have  given  clearly  defined  interference  fringes  in  the 
shock  layer  of  a  cylinder  immersed  in  a  shock-heated  flow  of  argon 
plasma.  Tlie  masking  effects  of  the  Intense  shock  layer  radiation  have 
been  elixrlnated  by  use  of  a  Kerr  cell  shutter  on  loan  through  the 
courtesy  of  the  nearby  Lockheed  Research  Laboratory.  The  optical  ex¬ 
periments  have  been  supplemented  by  a  series  of  preliminary-  thick  film 
heat  transfer  measurements  along  the  stagnation  line.  The  present  vork 
alms  to  measure  both  the  convective  and  radiative  heat  transfer  from  the 
shock  layer-boundary  la^/er  plasma,  with  independent  measurement  of  the 
electron  density  as  a  control.  As  such  it  will  represent  an  ex^^ension 
of  the  wor.r'  of  Rutowski  and  Bershader. 

Por  reference  purposes,  the  following  listing  of  output  of 
the  shock  tube  phase  cf  the  program  is  given: 

Technical  Supplements  to  Progress  Reports: 

H.  Wong  and  D.  Bershader,  "Interferometric  Study  of  Thermal  Equilibration 
of  a  Shock  Heated  Plasma,"  Supplement  2,  April  20,  196^^. 

R.  S.  deVoto,  "Transport  Properties  of  Partially  Ionized  Monatomic  Gases," 
being  Issued  as  a  supplement  to  this  report. 


Rutowski,  R.  and  Bershader,  D.,  Physics  of  Fluids,  568- 

577,  196^. 
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Teeioleai.  Society^  Presentatiocj 

H.  Ifeag,  «interfJro^tric  Stirfy  of  Shook- Induced  Reaction  Rates  In 

MWn,  AIAA  Aonxal  ^teeting,  ¥a?hingtoq,  D.C.,  July  1,  196^^. 


Journal  ftibllcatlons? 

The  work  of  both  ¥ong  and  deVbto  vill  appear  as  publicapions 
in  an  archive  Journal. 


Gi^duate  Thesis  Progr6iffisj 

&-ward  ^55!^  Ki.D.  (awned  June  1962^),  "Interferometric  Study  of 
Thermal  Equllibmtion  of  a  Shock  Heated  Plasma.” 

R.  S.  deVoto,^D.  (to  be  co^leted  January  1965),  "Transport  Properties 
of  Partially  Ionised  Monatomic  Gases.” 

The  following  thesis  programs  were  Independently  supported  but 

Bade  use  of  the  shock  tube  facility  and  dealt  with  the  closely  related 
subjects; 

Charles  ElUot-t,  Engineer’s  Degree  (completed  August  1964),  "Interfero- 

^  Inmiersed  in  a 

VmiamY.  Cole  E^lneer’s  Degree  (completed  August  1964),  ’^Measureme^t 

^  n  Shcck-Heated  Plasm  and  Comparison 
with  the  oaha  Equation.” 


All  in  ail,  +he  program  to  date  has  been  a  fruitful  one,  and 
it  will  be  coutinued  as  a  long-range  urogram  emphasizing  quantitative 
optical  nffiasurfimenta  of  dense,  shock-heated  plasmas. 


APPENDIX  A 


Spectral  Line  neversai  Method  for  Equilibi’iiim  Flows 

I,  Introduction 

The  sodium  line  reversal  method  of  flame  temperature 
measurement  was  introduced  by  Firy  in  1903*-^  Its  validity 
established  by  Kohn  In  1914^1  ana  engineering  measurements 
of  combustion  temperatures  well  underway  in  1933^^  Numerous 
references  describe  the  basic  theory  and  experimental  design 

T  ^  *1 

of  the  method,  some  of  tiie  best  are  references'^  ,  and  later 
1 H 

by  Mak  .  The  apparatus  developed  for  the  flame  temperature 
measurements  in  the  Stanford  MKD  generator  is  shown  in  Figure 
A-1.  A  tungsten  filament  lamp  shines  through  the  flame  con¬ 
taining  sodium  atoms,  whose  temperature  is  to  be  measured, 
and  then  passes  through  the  slit  of  the  spectrometer. 
Adjustment  of  the  lamp  orightness  tem.perature  to  the  tempera¬ 
ture  of  the  flame  will  cause  the  NaD  lines  (33 -3p  transition) 
to  match  with  the  lamp  background.  Precise  determi.nation  of 
the  calibrated  lamp  current -will  then  yield  tne  flame  tempera¬ 
ture  after  the  appropriate  corrections  for  brightness  tempera¬ 
ture,  wavelength,  and  optical  effects  have  bee.n  made. 

For  this  method  of  temperature  measurement  to  be  valid 
there  are  several  assumptions  which  must  be  made  about  the 
flame.  The  first  Is  that  the  flame  must  be  optically 
transparent  In  the  wavelength  region  of  interest  (589O-5696 
for  NaD  lines).  This  necessitates  the  use  of  small  amounts 
of  sodium  In  the  flame.  Secondly  the  flame  must  be  in 
thermodynamic  equilibrium  If  flame  temperatures  are  to  be 
measured.  Othei*  assumptions  are  listed  as  they  occur. 

II .  Sodium  Line  Reversal  Method:  T.heory 

* 

a .  Basic  Squat lens  and  Nomenc lature 


General  reference  for  this  section  is  18 


1)  The  radiation  density  in  energy  per  cm'  of  a 
black  body  at  temperature  T  ,  '’K  ,  having  wavelengths  bet'ween 
A  and  A  +  dA  can  be  formulated  by  Planck's  law 


SirhcA’^ 

--  hi  -1  ' 


exp 


kT 


exp 


c^A 
AT 


-5 


A 


(1) 


are  the  first  and  second  radiation  constan 


where  c^  and  c, 

c^  =  4.992  X  10~^^  erg  cm^  =  1.439  c:n°K. 


2}  The  radiation  intensity  in  ergs  per  sec 

(watts)  per  unit  solid  angle  normal  to  the  surface  of  the 
blackbody  of  area  A  is 

-5 


d  = 


2k  c^A 


(2) 


exp 


AT 


-  1 


3)  Wien's  law  is  a  very  good  approximation  for 

/  l8 

>5  in  the  visible  range  of  the  spectrum  so  that 


c, 

=  2Ac^  A‘^  exp  (-  j^) 


(3) 


4)  The  radiation  flux  over  a  solid  angle  O  is 


;hen 


(4) 

4a)  For  a  non-blackbcdy  v/ith  a  spectral  emissivity 
the  radiation  flux  is  reduced  to 

(5) 

5)  The  brightness  temperature  of  a  light  source 
(or  non-black  body)  at  a  wavelength  A  is  the  temperature 
at  which  it  is  necessary  to  .h^at  a  blackbody  so  that  its 
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brightness  matches  that  of  the  light  source  at  the  same 
wavelength.  The  brightness  temperature  can  be  deduced  from 
Equations  (3)  and  (4a)  above  as 


where  Is  the  spectral  emissivity  at  wavelength  A  cf 

the  eralttir.g  body  under  consideration,  T^.  Is  the  true 
temperature,  and  is  the  brightness  temperature. 

6)  Correction  of  light  source  brightness  temperature 
to  different  spectral  values  is  frequently  necessary  in  NaD 
line  reversal  methods.  This  is  accomplished  by  means  of  Eq. 

(6)  above.  Let  A  be  the  wavelength  of  the  spectrum  line 
used,  i.e.,  for  NaD  lines  X  =  O.589U.  ,  and  let  be 

the  calibration  wavelength  of  the  light  source.  Then 


_i - J.  .  in{c,) 


T  (A^) 
o  e 


Subtracting  gives 


1 


in(e  J 


7)  For  the  tungsten  filament  light  source  used, 
the  emissivity  must  be  known.  Larabee^^hs  developed  empirical 
equations  in  different  wavelength  regions  by  assuming  linear 
temperature  a.nd  wavelength  dependence.  For  A  =  O.450  to 
0.680U. 

€.  T,  »  0.4655  +  O.1558A  +  0.2675  -  lO’^T  -  0.7305  X  lO’^AT  (8) 
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8)  The  brightness  temperature  required  for  com¬ 
parison  In  the  NaD  line  reversal  method  is  that  of  the  aight 
source  at  its  image  position  in  the  flame.  It  is  thus  necessary 
to  correct  for  the  absorption  and  reflection  loss  of  light  by 
any  lenses,  windows,  or  filters  in  the  system  up  to  the  flame 
image  position.  (Such  losses  after  the  light  has  passed 
through  the  flame  merely  reduce  the  Intensity  of  bor.h  tne  NaD 
lines  and  the  background  light  source  radiation  simultaneously.) 
A  good  assumption  for  must  neutral  glasses  Is  that  the  spectral 

transmission  is  constant  except  in  the  ultraviolet  and  infared 

1  ft 

regions  where  it  drops  off  rapidly  .  If  I-^  is  the  transmis¬ 
sivity  of  the  glass,  the  emergent  flux  will  be 


emergent  ^A  ^A, incident  (9) 

From  this  relation,  substitution  for  0^  from  Planck's  law^ 
can  be  made  to  obtain  the  final  corrected  brightness  temperature 
Tg  of  the  source  at  the  image  position  in  the  flame  as 


Tb(a)  tyx) 


-f  s  +  a 

lens  windov/i 


+  a 


filters 


(10) 


wnere 


a 


(11) 


is  the  pyrometric  absorption  expressed 


b.  Theoretical  Basis  of  the  Sodium  Line  Reversal  Method 
The  basis  of  the  NaD  line  reversal  theory  is  developed 
below  first  by  the  energy  balance  method  for  clarity,  and  the 
theory  is  then  developed  from  the  equation  of  transfer  for  a 
better  understanding  of  the  physical  assumptions  invol'^ed. 

1.  Energy  Balance  method 

The  light  beam  from  the  light  source  at  a  brightness 
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temperature  Tg  passes  through  a  gaseous  body  at  temperature 
Tq  .  The  radiation  flux  is  then  specularly  isolated  by  the 
spectrometer,  received  by  a  photomultiplier  tube,  and  dis¬ 
played  for  measurement  by  a  suitable  Instrument.  The  radiation 
flux  thus  measured  is  the  sum  of  the  spectral  emission  of  the 
gas  Id  the  continuous  lamp  radiation  transmitted  by  the  gas. 
Thus, 


* 


By  assuming  thermal  equilibrium  and  applying  Kirchoff's  law 
=  a^)  ,  =  (l  -  a^)  the  above  becomes 

-  ^a)  -Tfi)  (^3) 

The  reversal  condition  is  that  the  total  radiative  flux  is 
unchanged  in  the  presence  of  the  gaseous  body,  when  this 
happens,  l.e.,  when  »  then  Tq  =  Tg  ,  or  the 

temperature  of  the  gas  is  equal  to  the  corrected  brightness 
temperature  of  the  light  source.  When  T^  is  greater  than 
Tg  >  is  larger  than  and  the  NaD  lines  appear  in 

emission.  When  T„  is  less  than  T^  ,  0-  .  .  is  less  than 
and  the  NaD  lines  appear  in  absorption.  When  T^  =  T.., 
the  lines  are  neither  in  absorption  nor  emission  but  would 
appear  to  be  a  continuum  on  a  low  resolutio.n  spectrometer. 
When  there  is  a  cool  boundary  layer  surrounding  a  central 
core  region  of  higher  temperature  the  lines  will  be  a 
combination  of  envnlssion  and  absorption  as  viewed  by  a  high 

# 

Many  isferences,  4,6,  for  example,  wrongly  omit  t.he  gas 
emissivlty  before  the  T^)  in  Equation  12.  Thl^ 

would  yield  the  result  that  tne  gas  emissivlty,  or 

absorption  coefficient  la  required  before  the  gas  tempera¬ 
ture  could  be  determined  from  the  NaD  line  reversal  method. 
This  conclusion  is  not  true,  as  t.he  follow l.ng  section  will 
demonstrate . 
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resolution  spectrometer.  The  observable  criteria  for  reversal 
is  then  that  the  Integral  of  the  observed  flux  over  the  wave¬ 
length  region  of  the  lines  be  zero.  Reference  16  discusses 
this  effect  in  detail. 

2.  Solution  of  the  Equation  of  Transfer^^ 

Assumptions ; 

a)  The  flow  is  one  dimensional  and  planar. 

b)  The  gas  is  in  local  thermodynamical  equilibrium. 

c)  The  gas  is  optically  transparent. 

d'  The  presence  of  small  amounts  of  sodium  does  not 
affect  the  properties  of  the  test  gas. 

e)  The  electronic  excitation  of  the  sodium  atoms 
is  in  equilibrium  with  the  test  gas. 

The  equation  of  transfer  under  these  assumptions  is 


—  =  k,  p(B,  -  Ij 


(14) 


where  is  the  mass  absorption  coefficient  of  the  gas, 

is  the  Planck  function*,  and  p  is  the  density -of  the  gas. 
(a  Vkj^p  is  the  radiation  mean  free  path  of  the  gas). 
Solving  Eq.  (l4). 


(l4a) 


where  n  =  -  B^  and  is  the  optical  depth  of  the  gas. 


3,  =  /  h.  p  d  X 


From  Eq .  ( 14a ) , 


h  -  h  '  =  ® 


(14b) 


is  derivable  from  Eq.  1  by  use  of  the  relations  Av  = 

ana  -(o/47r)u^  ;  -  2h>.--c,,  [  esp(h,/h  )  -  IJ  , 


c 


Applying  the  boundary  condition  that  (lamp  intensity) 


the  above  becomes 


-s. 


«  )  +  h,B  ® 


-S- 


(l4o) 


Integrating  Eq.  (I4c)  over  the  wavelength  region  of  the  NaD 
line  gives 


-s. 


Lr.  ■  A.  -  -  e  '')dX  +  /  I 


-3. 


NaD 


NaD 


NaD 


dX  (Ikd) 


A  typical  NaD  llneat  the  reversal  condition  is  shown  below 
as  seen  by  a  spectroscope  of  high  resolving  power  through 
a  gas  with  cool  edges. 


NaD  Line  at  Reversal  Condition.  Tne  reversal 
condition  is  that  +  A^  =  A^  . 

Then,  assuming  symmetrical  waveieiigth  distribution 

Ld  °  =  2w  g  (14e) 

Taking  independent  of  X  over  the  shcrt  region  of 

integration  gives 


^X,B  *  ^  ^  +  ^X  B  ®  ^ 

NaD  j^aD 
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or 


■A,B 


-  ej 


-S. 


2vi 


/  e 
NaD 


dX 


(I4g) 


rearranging  gives 


S  “  s 

2w  NaD  ^  NaD 

Thus  the  condition  for  reversal  gives 


^X.B  =  \  Tg  =  Tg  (141) 

independent  of  the  optical  depth  Sj^  . 
above  analysis  contains  the  additional  assumptions 

a)  =  constant  with  respect  to  wavelength  over 
the  NaD  line  wavelength  region. 

b)  X  -  X  =  X  -  X.  ;  l.e.,  line  symmetry. 

c)  “X  B  ~  constant  with  respect  to  wavelength  over 
the.  NaD  line. 

3.  Cool  Boundary  Layer  Effects 

In  actual  flow  in  an  MHD  generator  there  is  a  boundary 

layer  growth  along  the  walls.  Any  sodium  vapor  trapped  in  this 

cool  layer  will  absorb  radiation  from  the  emitted  radiation  of 

the  central  core  region  of  hot  gas.  The  quantities 

0^  measured  will  then  be  some  mean  quantities  along  the  optical 

path  length  in  the  test  section.  The  gas  temperature  measured  is 

thus  a  mean  temperature  T„  which  is  less  than  the  true  free- 

G,  mi 

stream  value  T„  .  T_  will  depend  on  the  boundary  conditions, 

G,  t  G, m 

thermal  gradients  present  in  the  cool  boundary  layer,  and  the 
number  density  of  sodium  atoms  present. 

Consj.der*  a  flow  consisting  of  three  separate  temperature 
regions,  designated  1,  2  and  3.  The  solution  of  the  equation 
oh  transfer  through  o;.  'f  these  layers  is  given  by  Eq.  (l^.c) 
as 


which  is 
The 

that 
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^  -  e 


(15a) 


where  Ig  is  the  radiation  intensity  Incident  on  layer  1 
(corresponding  to  lamp  intensity)  am  is  the  intensity 

leaving  layer  1  and  passing  on  to  layer  2.  Likewise 


)  +  B^d 


(l5o) 


where  is  the  intensity  leaving  layer  2  and  passing  on  to 

layer  3. 


h: 


ot 


B^d-e 
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B^d-e 


-s, 


)e 


-Sn 


B^(l-( 


-S'. 


+  I 


"2  + 
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S3) 


(15c) 


By  taking  layers  1  and  3  to  be  boundary  layers  with  identical 
temperatures  and  optical  depths,  and  assuming  and  s,, 

snuill,  Eq.  15c  can  be  solved  for  the  core  condition  is 


(1  -  e  “2)  (1  -  3^) 

The  observed  mean  value  is  B^  ,  where  for  reversal 

(jT  ^  m 

tions  =  ig  .  Thus  the  relation  between 

temperature  and  the  mean  temperature  observed  is 


d5d) 


condl- 
the  core 


-3, 


(1  -  asp)  - 


ipid 


(I5e) 


For  order  of  magnitude  estimates  which  do  not  Involve  the 
optical  depths  ox’  ,  the  following  equation  is  convenient^^ 
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Ill .  Sodium  Line  Reversal  Method;  Experimental 

a.  Apparatus 

1.  Light  Source 

The  light  source  used  In  the  present  experimental 
facility  (see  Figure  1)  Is  a  calibrated  General  Electric 
tungsten  filament  lamp,  model  number  i8A(T10)1.  The 
calibration  cavers  the  true  temperature  range  1‘rom  2600  K  to 
3200^K.  The  lamp  power  supply  is  a  motor  driven  three  phase 
varlac  and  rectifier  which  provides  up  to  25  amperes  of- 
rectified  direct  current.  Two  27  ampere  rheostats  are 
connected  In  parallel  to  provide  fine  current  control.  The 
current  to  the  lamp  is  determineo  by  precisely  measuring  the 
voltage  drop  across  an  accurate  25  ampere  shunt.  Due  to  the 
thermal  lag  of  the  lamp  filament,  the  rectification  ripple  is 
not  observable  when  the  lamp  radiation  Is  measured  ;vith  the 
spectrometer  appratus  utilized.  The  power  supply-lamp  system 
has  pr‘''-on  to  be  extremely  stable  and  thus  reliable  for  the 
spectral  line  reversal  method. 

2.  Spectrometer  System 

A  Jarrell-Ash  half-meter  Ebert  scanning  spectrometer 
has  been  temporarily  borrowe'd  for  determination  of  spectral 
intensities.  This  instrument  has  proved  to  be  highly  success¬ 
ful  for  the  spectral  line  reversal  measurements.  The  Jarrell- 
Ash  power  supply  for  the  photomultiplier  and  the  photomultiplier 
output  amplifier  for  the  above  spectrometer  are  combined  Into 
one  unit,  the  output  Is  recorded  on  a  Textronlcs  Type  5^1 
oscilloscope  with  maximum  sensitivity  of  1  mv/centlraeter . 

The  spectral  line  profile  Is  obtained  by  photographing  the 
sweep  of  the  oscilloscope  ever  the  spectral  line  wavelength 
region  by  a  time  exposure  with  the  Polaroid  attachment  provided 
with  the  oscilloscope.  The  spectral  line  can  also  be  observed 
directly  or.  the  oscilloscope. 

3.  Lenses 

The  lenses  utilized  for  the  present  apparatus  are 
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two  glass  7  centimeter  diameter^  20  centimeter  focal  length 
identical  lenses.  Their  arrangement  In  the  optical  design 
is  dictated  by  the  Internal  optical  arrangement  of  the 
Jarrell-Ash  spectrometer.  The  optical  design  is  shown  in 
Figure  1. 

b)  Experimental  Observations 

Typical  examples  of  the  NaD  lines  as  seen  in  the  small 
propane  flame  are  shown  on  Figures  i,  4,  5,  and  6.  The 
temperature  of  the  flame  is  determined  to  be  2135 °K  as  seen 
from  Figure  6^  bottom.  The  temperature  measurements  have 
been  determined  to  an  observable  experimental  accuracy  of 
10°K.  By  this  it  is  meant  that  two  different  output  spectral 
line  shapes  corresponding  to  a  difference  of  lO^K  in  flame 
temperature  are  easily  distinguished  and  accurately  measured. 
Measurements  on  the  MHD  channel  at  the  plenum  chamber  exit, 
exhausting  to  the  atmosphere,  have  shown  good  results  with 
the  existing  apparatus.  Preliminary  results  show  excellent 
line  shapes  and  reversal  lines  with  a  sodium  seed  concentra¬ 
tion  of  less  than  0.001^.  The  temperatures  determined  by 
the  spectral  line  reversal  method  agree  ;vell  with  the 
theoretical  calculations  for  the  particular  flow  rates 
chosen. 
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FIGURE  A 


FIGURE  A-3 


2  5890  and  5896  A  NaD  lines  in  emission. 

Light  source  off.  10  millivolt  per 
division,  0.2  seconds  per  division. 


5890  and  5896  A  HaD  lines  in  absorption. 
15*80  amperes  lamp  current. 
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FIGURE  A“4  NaD  lines  near  reversal  condition  but 
slightly  in  emission.  12.65  amperes 
lamp  current,  corresponding  to  a  gas 
temperature  of  2125‘*K. 


FIGURE  A--5  Top:  NaD  lines  near  reversal  conditions. 

12.50  amperes,  2110®K. 

Bottom:  NaD  lines  at  reversal  conditions 
12.78  amperes.  Gas  temperature 
determined  to  be  2135‘’K. 
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n.  ELECTROCHEMICAL  CONVERSION 


A.  The  Effects  of  Free-Stream  Turbulence  on  the  Limiting  Current  Density 
at  the  Front  Stagnation  Point  of  a  Translating  Cylinder  Electrode 

Introduction 

The  overall  rate  of  conversion  in  an  electrochemical  reaction  may  be 
controlled  by  either  of  two  processes*  the  rate  of  conversion  at  the  electrode 
surface,  or  the  diffusion  of  the  reactants  through  the  boundary  layer.  For 
the  systems  v/hich  are  chemically  slow-reacting,  the  diffusion  process  is 
relatively  fast  and  the  electrode  surface  concentration  of  reactants  is  equal 
to  the  bulk  concentration.  The  bulk  concentration  would  be  a  known  parameter 
of  the  system. 

Controlled  hydrodynamic al  techniques  are  used  to  study  the  very  fast 
electrochemical  reactions,  where  the  process  maybe,  at  first  diffusion 
controlled.  If  the  rate  of  diffusion  may  be  systematically  controlled,  it  may 
be  possible  to  approach  a  regime  where  the  overall  rate  is  again  governed  by 
the  kinetics.  If,  in  addition,  it  is  possible  to  express  the  diffusional  process 
from  a  theoretical  standpoint,  it  would  be  possible  to  study  the  kinetics 
during  the  intermediate  regime  since  the  surface  concentration  would  be 
known  from  theoretical  considerations.  It  would  be  impossible  to  say  anything 
about  the  surface  concentration  in  a  completely  diffusional-controlled  system. 

The  rotating  disk  electrode  and  the  translating  cylinder  electrode 
have  been  used,  successfully,  to  study  the  faster  electrochemical  reactions. 
Both  of  these  devices  offer  an  electrode  surface  of  uniform  concentration;  if, 
in  addition,  the  entire  surface  is  maintained  at  a  constant  potential.  The 
laminar  boundary  layer  theory  for  the  rotating  disk  was  first  presented  by 
Levich  (4).  The  theory  for  the  translating  cylinder  and  its  possible  use  as 
an  electrochemical  tool  was  first  suggested  by  Bopp  and  Mason  (Z). 

The  interest  here  lies  in  the  effect  of  the  free  stream  turbulence 
gene  rated  by  the  cylinder  vipon  the  laminar  boundary  layer  at  the  front  stagnation 
region.  It  must  be  demonstrated  that  the  laminar  boundary  layer  theory  is, 
indeed,  applicable,  and  to  within  what  experimental  limitations  it  may  be 
applied. 
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Experimental  Apparatus  and  Procedure 


With  the  exception  of  a  few  modifications,  the  experimental 
arrangement  which  was  used  here  corresponds  to  that  used  in  the  studies  by 
Bopp  and  Mason  (2).  The  tank  was  a  standard  cylindrical  pyrex  jar,  25.  4  cm. 
high  and  25.  7  cm.  O.  D.  ,  closed  at  the  top  by  a  1.  9  cm.  thick  Incite  flange  . 
The  tank  rested  on  an  aluminum  plate  which  contained  four  vertical  studs. 

The  studs  rose  up  through  the  top  flange  with  the  flange  then  being  secured 
to  the  tank  by  four  nuts  (See  Fig.  2,  Ref.  2  ).  A  Viton  A  ring  gasket  offered 
an  effective  seal  between  the  Incite  flange  and  the  top  of  the  pyrex  tank.  Four 
vertical  baffl''"  were  placed  inside  the  tank  and  connected  to  the  top  flange 
by  nylon  screws.  Any  number  of  the  baffles  could  then  be  removed  if  it  was  so 
desired.  The  bearing  housing  and  shaft  arrangement  is  shown  schematically 
in  Fig.  7  of  Bopp  and  Mason  (2).  The  shaft  rotated  on  two  bearings,  3.  8  cm. 
apart.  The  bearings  were  sealed  from  the  solution  in  the  tank  by  a  Teflon- 
lined  Viton-A  annular  cylinder,  0.  ?5  inch  I.  D.  by  1.  0  inch  O.  D.  and 
0.  5  Inch  thick.  The  entire  essembly  was  attached  to  the  top  flange  by  six 
screws  at  the  base  of  the  bearing  housing.  The  shaft  was  coated  with  an  epoxy 
resin  to  reduce  the  damage  due  to  chemical  attack.  A  Heller  type  GT21 
laboratory  mixer  and  controller  was  used  to  rotating  the  shaft.  By  varying 
the  pvilly  ratios  between  the  motor  and  the  shaft,  rotational  speeds  of  from 
2  rpm  to  170  rpm  could  be  maintained  and  controlled  quite  successfully.  The 
auxiliary  electrode,  a  3.  57  cm. by  4.  45  cm.  piece  of  platinised  platinum,  was 
placed  on  the  base  of  the  tank.  A  60  cm.  length  of  25  ml.  platinum  wire  was 
<;pot-welded  to  the  auxiliary  electrode  and  drawn  up  along  one  of  the  baffles 
and  out  the  fuel  port.  The  lead  from  the  st  electrode  in  the  cylinder  was 
drawn  up  through  the  hollow  shaft  into  a  mercury  dish,  which  rotated  on  the 
shaft.  The  circuit  was  then  completed  by  dipping  another  external  lead  into 
the  dish.  A  lugin-type  calomel  electrode  was  used  as  a  reference  electrode. 
The  capillary  tube  extended  down  into  ne  supporting  electrolyte  through  the 


top  flange.  A  gas-bubbler  was  also  supplied  through  the  top  flange.  Nitrogen 
gas  was  continually  swept  through  the  solution  to  prevent  the  build-up  of 
dissolved  oxygen  in  the  supporting  electrolyte.  Three  Teflon  Swagelock 
fittings  were  used  in  conjunction  with  the  three  ports  in  the  top  flange.  A 
model  855B  power  supply  from  the  Harrison  Laboratories  was  used  to  control 
the  potential  between  the  test  and  the  reference  electrodes.  A  J-Omega  type 
21A  D.  C.  volt-Qhm  meter  was  used  to  monitor  the  applied  potential.  The 
current  was  recorded  on  a  Moseley  Autograph  X-Y  Plotter.  The  rotational 
speeds  were  measured,  to  within  an  accuracy  of  about  2  °/o,  by  a  hand  stop 
watch. 

The  electrode  surfaces  were  cleaned  in  a  strong  solution  of  KOH  for 

periods  up  to  one-half  hour.  They  were  then  used  and  placed  in  a  50  °/o 

aqua  regia  solution  for  fifteen  minutes.  Finally  the  electrodes  were 

cathodized  for  one  to  two  hours  in  0.  IN  H,SO^  at  100  before  being  placed 

2  4  cm^- 

into  the  test  solution. 

Procedure ; 

The  voltage  was  first  adjusted,  either  cathodically  or  anodically, 
until  the  system  was  operating  on  the  current  plateau  (Fig.  1).  Once  a 
rotational  speed  was  selected,  three  minutes  were  allowed  before  the  X-Y 
Plotter  was  engaged.  After  each  run,  five  minutes  elapsed  before  the  cylinder 
was  again  rotated.  Any  turbulence  which  had  been  previously  generated  was, 
by  then,  dampened  out. 


Discussion 


Electrochemical 

Consider  the  electrochemical  reaction  occurring  between  species 
R  and  Ox; 

R  ^  Ox  +  ne" 

where  n  electrons  are  transfered  per  mole  of  reactant.  One  may  write  the 
rates  of  reaction  for  the  forward  and  reverse  steps,  assuming  for  example 
1st  order,  as: 


where 

r,  ,  r,  -  rate  of  reverse  and  forward  reaction 
b  t 


is) 

Ox 


surface  concentration  of  R  and  Ox 


j  -  rate  constant 
o 


oc  -  transfer  coefficient 

F  -  Faradays  constant  96,500  - 

equivalent 

n  -  number  of  equivalents  per  mole 
R  -  gas  constant 
T  -  temperature 
V  -  overvoltage 


At  the  equilibrium  rest  potential,  ^  =  0.  The  sign  of  r)  determines  whether 
one  is  operating  cathodically  or  anodically.  As  the  overvoltage  is  increased 
anodically,  the  rate  of  the  forward  reaction  increases  with  respect  to  the 
reverse  reaction  (Eq.  1).  The  current  density  measured  at  the  electrode 
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will  then  increase  from  zero  and  reach  a  plateau.  The  plateau  will  occur 
when  (r^  -  r^^)  is  sufficiently  large  that  the  system  becomes  diffusion 
controlled.  The  ratio  of  the  forward  to  reverse  rates  of  reaction  may  be 
formed  directly  from  Eq,  1: 


nnF 
>•  ‘  ' 
RT 

m 


In  the  diffusion  controlled  regime,  the  value  oi  C'  is  at  least  tv/o  orders 

iv 

of  magnitude  less  than  the  bulk  concentration.  If  the  two  species  have 
-  •  (s) 

Schmi'^t  numbers  of  the  same  oraer  of  magnitude,  C'  '  will  be  about  double 

Ox 

the  bulk  value.  Thus,  if  the  solutions  are  equimolar  in  the  bulk,  for  n  =  1 


10*’ C' 


:”3»  V 


=  10“’ e-^*^ 

(crt\ 

,  The  surface  concentration  of  R  may  be  less  than  10’ ’C'  '  but  an 

IV 

anodic  overvoltage  of  0.  300  -  0.  400  volts  (77  is  negative)  is  usually 
s^lfficient  to  reach  the  plateau.  Figure  1  shows  a  polarization  curve 
obtained  for  the  oxidation  of  the  ferrocyanide  ion  to  the  ferric yanide  ion  in 
1  M  KCl. 


Hydrodynamic  s 

The  co-ordinate  system  which  is  used  for  two  dimensional  boundary- 
layer  flow  past  a  cylinder  is  chosen  such  that  the  positive  y-axis  extends  out 
away  from  the  cylinder  in  the  radial  direction  and  the  positive  x-axis  increases 
clockwise  in  the  circumferential  direction  from  the  front  stagnation  point.  Of 
interest  here  is  the  asymptotic  case  of  Re  -  00,  Sc  -»  00.  The  flow  and 
continuity  equations  may  be  written  in  dimensionless  form  with  all  velocities 
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relative  to  and  all  lengths  relative  to  the  cylinder  diameter,  d: 


du  .  3u 

u  ““  T  V  ~~ 

dx  dy 


du  ^ 

dx  dy 


JL  Phi 

;ie  >y2 


-  U(x) 


dU(x) 


For  flow  past  the  front  stagnation  point  of  a  cylinder  U(x)  =  x  (4).  Letting 
u  =  xf  {p),  withTj  =  '/^  y,  Eq.  3a  may  be  transformed  to 

f"  +  ff"  +  (1  -f  2)  =  0  (3b) 

Equation  3b  has  been  solved  numerically  and  f(t^)  is  tabulated  in  Schlichting  (4). 


Defining  a  dimensionless  concentration,  0= — 


;re 


{ s ) 

is  bulk  and  C  is  the  surface  concentration,  one  may  write  the  mass 
transfer  equations  as: 


bd  .  „  39 

u  T  V 

dx  bv 


1  b’^s 
Sc  dfi* 


where 


Re  V 


Boundary  conditions 

T7  =  0  0=1 

p  -*  <X)  6  —  0 

For  Sc  -♦  00,  the  mass  boundary  layer  is  much  thinner  than  the  momentum 
boundary  layer.  One  may,  then,  expand  u  linearly  in  rj.  Thus: 


p  =  fiMrj 


p=  0 


6l 


and 

-I  i  [H  >7’  (5) 

One  may  eliminate  Sc  fvom  Eq.  4  by  istting  z  =  Sc  Under  the  above 

transformations,  one  "htains  the  new  equation  (1); 


a^9 


+ 


=  0 
34 


Boundary  Conditions 
C  =  0  0=1 

^  OO  0  -«  0 


(7b) 
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Solving  Eq.  7b  under  the  imposed  boundary  conditions 


If  e  dt 


■<f) 


The  diffusional  flux  of  the  i  component  to  the  surface  may  now  be  easily 
calculated. 


=  -D.  p- 

im  dy  - 
'  y  =  0 


with 


S/x)  =  f"(0)  U(x)  =  1.2326X 


Thus 


D.  rc<“».C<^| 
im\  A 


[u.660Re^/®Sc^/®] 


The  diffusion -controlled  current  depiv,  .  J  by  the  plateau  of  Fig.  1  corresponds 
to  the  condition  (c<7.  0.  The  current  is  related  simply  to  the  concentration 
gradient  at  the  electrode  by: 


I  =  -nFD. 


im  dy 


y  =  0 


where 


n  =  charge  of  i  component 


F  =  96, 500 


coulombs 

equivalent 
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For  the  translating  electrode,  Eq.  9  may  be  put  into  a  more  convenient  form. 
Let: 

R  V 

-  8  ^ 


where 

O  =  rotational  speed  (rpm) 

R  =  shaft  radius 

s 

Thus 

I.  *  (4. X  10*  )ni;H  D.  (11 

^  im  y  d  /  i  ' 

Figs.  4,  5,  and  6  are  graphs  indicating  the  linearity  of  I  with  for  the 
ferricyanide-ferrocyanide  system. 

The  theory  presented  above  ia  for  free  stream  flow  past  a 
stationary  cylinder.  It  is  more  practical,  however,  to  rotate  the  cylinder 
in  a  stationary  bed  of  fluid.  The  equations  presented  above  are  applicable 
to  both  'low  systems  if  the  boundary  conditions  are  the  same.  Translating 
the  cylinder  through  a  stationary  fluid  does,  however,  create  a  :,ertain 
degree  of  free  stream  turbulence.  The  boundar/  conditions  at  the  outer 
edge  of  the  boundary  layer  will  not,  therefore,  be  identical  to  those  used 
in  the  formulation  of  Eq.  11.  The  use  of  laminar  boundary  layer  theory 
to  describe  the  translating  cylinder  must  be  justified  experimentally.  It 
is  clear  at  the  outset  that  there  will  be  certain  limitations.  The  size  of 
the  cylinder  will  determine  the  intensity  of  the  free  stream  turbulence, 
and  if  deviations  from  the  laminar  theory  occur,  they  will  be  predominant 
with  the  larger  cylinders. 


Experimental  Results 


The  potassium  ferricyanide-ferrocyanide  system  has  been  selected 
for  this  work  since  it  has  been  studied  by  other  investigators  and  has  been 
shown  to  be,  electrochemically,  a  relatively  rapid  reaction.  Of  initial 
interest  are  the  values  of  the  diffusion  coefficients  for  the  reduced  and 
oxidized  species.  These  values  have  been  calculated  by  other  workers  (3,  6) 
using  various  techniques  and  supporting  electrolytes.  The  values  which  have 
been  obtained  here  are  given  ir  Table  III.  The  values  have  been  obtained 
from  a  least  squares  analysis  >f  the  slopes  of  the  graphs  shown  Figs.  4,  5, 
and  6  in  accordance  with  the  laminar  boundary  layer  theory  presented  in  the 
previous  section.  These  grap'us  are  bound  in  a  region  of  Reynolds  number 
from  ZOO  to  8,000  and  are  well  vithin  the  limits  of  the  laminar  boundary 
layer  theory. 

The  tank,  as  mentior  ed  in  the  previous  section,  has  been  equipped  with  four 
equally  spaced  baffles  around  its  perimeter.  The  baffles  offer  resistance 
to  tlie  main  body  of  the  fluid  and  reduce  the  overall  rotation  of  the  bulk. 
Observations  of  small  particles  in  the  fluid  indicated  that  at  least  until  the 
point  where  free  stream  turbulence  began  to  control  the  transport  process 
the  bulk  rotation  was  negligible.  If  the  fluid  were  rotating  in  the  bulk  the 
apparent  Reynolds  number  would  be  greater  than  the  actual  value  and  in 
curves  similar  to  Figs.  7,8,  9,  and  10  the  values  of  the  current  density 
would  tend  to  fall  below  the  line  given  by  the  laminar  theory.  This  was, 
however,  never  observed. 

Figure  3,  which  is  a  graph  of  the  frequency  of  turbulence  as  a 
function  of  rotational  speeds  indicates  that  turbulence  may  be  present  in 
the  tank  at  relatively  low  Reynolds  numbers.  The  frequency  increases 
rapidly  at  the  lower  rotational  speeds  and  then  tends  to  level  off  with 
increase  speed  of  rotation.  The  graph  shows  points  taken  with  all  of  the 
cylinders,  and  with  the  exception  of  the  large  electrode  on  the  1"  cylinder, 
the  frequencies  are  of  the  same  order  of  magnitude  independent  of  Reynolds 
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number,  and  a  function  only  of  the  rotational  speed.  The  1”  cylinder 
gave  lower  frequencies  for  a  given  speed.  The  amplitude  of  the  fluctuation 
is,  however,  a  marked  function  of  the  cylinder  size,  decreasing  with 
decreasing  diameter  of  the  cylinder.  This  is,  indeed,  what  one  might 
intuitively  expect  if  it  is  the  cylinder  itself  which  is  creating  the  free 
stream  disturbances  which  are  in  turn  causing  the  current  fluctuations. 

The  variation  of  current  amplitude  with  c /Under  size  is  given  in  Table  IV. 
The  current  fluctuations  which  were  recorded  are  direct  measures  of  the 
fluctuations  in  the  concentration  gradient  through  the  mass  boundary  layer. 
Since  C<oc) 

is  a  constant  and  0,  the  fluctuations  indicate  that  the 

boundary  layer  thickness  is  oscillating  due  to  the  influence  of  the  free 
stream  turbulence. 

As  would  be  expected,  the  size  of  the  cylinder  is  a  key  factor  in 
determining  the  turbulent  characteristics  of  the  system.  Figures  7  through 
16  are  graphs  of  the  current  and  the  square  root  of  the  rotational  speed.  The 
currents  have  been  taken  as  the  mean  of  the  fluctuations.  In  Figure  12, 
where  there  is  no  apparent  curvature,  the  linear  portion  extends  from  a 
frequency  of  0.  08  c.  p.  s.  to  about  3,  00c.  p.  s.  Since  the  line  extends  to 
such  small  frequencies,  where  the  laminar  theory  would  apply,  if  anN-where 
at  all,  it  can  be  argued  that  the  mean  of  the  fluctuation  is  given  by  the 
laminar  theory  up  to  the  point  of  curvature.  It  is  upon  this  basis  that  the 
linear  portion  has  been  used  to  determine  diffusion  coefficients. 

Equation  11  indicates  that  from  the  standpoint  of  the  boundary- 
layer  theory  it  would  be  advantageous  to  make  the  size  of  the  cylinder  as 
small  as  feasible.  Thus,  although  for  a  given  rotational  speed  the  Reynolds 
number  has  decreased,  the  concentration  gradient  has  increased.  As  will 
be  indicated  subsequently,  however,  there  is  a  practical  limitation  to  the 
cylinder  size  as  governed  by  the  size  of  the  electrode.  Figures  7,  8,  9, 
and  10  show  that  the  rotational  speed  at  which  curvature  begins  increases 
as  the  cylinder  diameter  decreases,  with  no  apparent  curvature  present 
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in  Figures  11,  12,  and  13.  The  increase  in  Nusselt  number  must  be 
attributed  to  free  stream  turbulence  and  not  premature  tripping  of  the 
momentum  boundary  layer.  Indeed,  curvature  has  been  observed  in 
the  region  of  Reynolds  numbers  from  3,  000  -  10,  000  where  the  boundary 
layer  should  be  quite  stable  with  respect  to  external  effects. 

It  can  also  be  shown  that  the  width  of  the  electrode  surface 

governs  the  onset  of  curvature.  Figures  7,  8,  and  9  depict  the  current- 

speed  behavior  in  the  same  tank  with  the  same  size  cylinder,  but  with 

different  size  electrodes.  When  one  av'erages  the  Nusselt  number  over 

a  larger  surface  area,  the  curvature  bct-ins  at  a  lov'jr  speed.  By 

o  o 

reducing  the  angle  of  the  electrode  from  8  to  less  than  1  the  speed  at 
which  curvature  begins  may  be  increased  by  45°/o.  This  effect  is 
brought  about  even  more  vividly  in  Figure  14  where  the  current-speed 
graph  is  given  for  the  0.  125"  cylinder.  Although  this  cylinder  was  small 
enough  to  cause  small  amplitude  fluctuations,  as  given  in  Table  IV,  the 
angle  of  the  electrode  was  14°,  large  enough  to  cause  curvature  at  a 
Reynolds  number  of  3,  000.  The  average  Nusselt  number  will  vary 
linearly  with  the  square  root  of  the  Reynolds  number,  and  independently 
of  the  size  of  the  electrode,  if  the  velocity  in  the  potential  flow  regime 
varies  linearly  with  x.  Since  curvature  appears  with  the  small  cylinder 
at  an  angle  of  14°,  the  free  stream  turbulence  introduces  higher  order 
terms  in  x  much  sooner  than  the  laminar  theory  would  predict.  For  the 
boundary  layer  theory  an  angle  of  20°  should  still  suffice  for  the  stagnation 
region.  Figures  7,  8,  and  9  with  the  I”  cylinder  indicates  that  the  higher 
order  terms  come  into  play  at  smaller  values  of  the  x-coordinate  when 
the  scale  of  turbulence  increases. 

It  is  also  of  interest  to  note  that  in  Figures  7  and  9  after  the 
transition  from  the  initial  straight  line  occurs,  the  graph  assumes  a  new 
line  with  a  larger  slopei  It  is  possible  to  explain  this  behavior  if  a  quasi¬ 
steady  state  mechanism  is  assumed  for  the  mass  and  momentum  transport 
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processes  in  the  presence  of  turbulence.  In  a  quasi-steady  state  sytem 
the  velocity  and  concentration  profiles  are  assumed  to  satisfy  the  time- 
independ  transfer  equations  (Eq.  3)  at  any  instance  of  time.  Under  the 
framework  of  the  boundary  layer  transformation  the  spacially  dependenc 
Nusselt  number  may  be  expressed,  for  any  geometry  or  potential  velocity 
field,  as 

Nu(x)  =  b(x)  Re  Sc  (12) 

This  is  the  asymptotic  expression  for  Sc  -*  oo.  For  stagnation  regime 
flows  b{x)  =  constant.  In  Equation  12  the  Nu(x)  has  been  time  averaged, 
and  the  average  value  is  assumed  to  be  time  independent.  Observations 
of  the  ciirrent  fluctuations  in  Figure  2,  indicate  that  the  mean  current 
is  independent  of  time,  supporting  this  assumption.  As  the  free  stream 
velocity  alters  its  spacially  dependent  character  with  changes  in  the 
intensity  of  turbulence,  b(x)  becomes  a  function  of  rotational  sp.::ed.  As 
mentioned  above,  the  frequency  of  turbulence  appeared  to  level  off  after 
its  initial  rapid  rate  increase.  If  this  may  be  taken  as  also  a  region  where 
the  turbulent  effects  level  off,  then  b(x)  would  become  only  a  function  of 
X  again.  Thus 

f*o  _  f  *o 

*0  Nu(x)dx  =  Nu  =  Re*'^®  Sc^'"^  4  b(x)dx 

^o 

=  Rc'/>Sci^3a(xJ. 

In  Figs.  7,8,  and  9  the  size  of  the  electrodes  were  constant,  thus 
a(Xo)  would  be  a  parameter  of  the  system.  The  system  would  then  become 
linear  again  with  the  square  root  of  Reynolds  number.  This  argument 
can  only  be  substintiated  more  conclusively  after  further  studies  on  the 
effects  of  the  electrode  size  on  the  average  value  of  the  Nusselt  number. 
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At  the  present  there  is  no  clear  explanation  for  the  difference  in 
the  observed  points  of  curvature  anodically  and  cathodically .  The  solutions 
were  equimolar  in  ferri-and  ferro-constituents  and  diffusion  coefficients. 
This  result  was  reproducible  in  a  series  of  runs.  It  might  be  possible  that 
an  explanation  lies  in  the  elucidation  of  the  effects  of  transporting  a  charged 
species  through  an  oscillating  velocity  field.  The  curvature  cannot  be 
attributed  entirely,  however,  to  the  charge  of  the  species  being  transported 
to  the  surface  of  the  electrode  since  curvature  was  also  obtained  with  the 
hydrogen  gas  ionization  reaction  as  shown  in  Figs.  15  and  16.  Here  the 
species  of  importance  is  the  neutral  H.,  molecule.  Also, decreasing  the 
size  of  the  electrode  increased  the  rotational  speed  at  which  the  onset  of 
curvature  was  noted  for  the  hydrogen  system.  This  is  shown  in  Fig.  16. 
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Conclusions 


Although  fte  froquoncy  of  turbulence  did  not  appear  to  be  a  function 
o  »  cylinder  size,  the  amplitude  did  decrease  as  the  cylinder  diameter 
ecreased.  The  1"  and  3/4”  cylinders  gave  curvature  whereas  the  1/2” 

and  1/4”  cylinder,  obeyed  the  laminar  theory  in  the  ra,.ge  of  rotational 
speeds  that  were  investigated  (2-170  rum) 

^  '  wu  rpm).  The  curvature  in  the  1/8" 

cywer  is  attributed  to  the  large  angle  ,14°,  which  the  electrode  subtends 
»  a  ront  face  of  the  cylinder.  Changing  the  angle  from  8°  to  1°  i,.  the 
cylinder  increased  the  speed  at  which  curvatu  e  commenced  by  45  °/o 
t  Should  be  observed,  though,  that  although  curvature  in  the  1/8”  cylinder 
occurred  at  a  Reynolds  number  of  about  one-third  that  of  the  1”  cylinder 
e  current  corresponding  to  the  point  of  curvature  was  about  three  times 
e»  erge  for  the  1/8"  relative  to  the  1”  cylinder.  1,  would,  of  course,  be 
^ssible  to  use  a  cylinde  i -electrode  arrangement  that  gave  curvature  if 
controlled  regime  remained  bounded  below  the  region  of 
^urvature.  It  still  must  be  verified  experimentally  that  when  one  approaches 
metics -controlled  regime,  that  the  presence  of  a  non-vanishing  surface 
concentration  does  not  initiate  deviations  from  the  laminar  theory  A.  was 
noted  above,  the  point  of  curvature  was  different  anodically  than  ca.hodica.ly 
-ggcsting  that  Uie  charge  of  the  species  plays  some  role  in  the  deviation 

rom  inearity.  I.  should  also  be  possible,  after  further  experimentation  to 

correlate  the  data  and  obtain  an  empirical  fit  for  the  second  r 

second  linear  portion  of 

the  graphs  in  Figures  7  and  9.  The  correlation  would  involve  varying  the 
aue  of  the  cylinder,  eleCrode.and  possibly  «,e  ratio  of  the  shaft  radius  to 

nh  radius.  An  empirical  relation.if  so  obtained,  would  greatly  increase 

tne  region  of  applicability  of  the  tran-?l-if t  j 

y  translating  cylinder  electrode. 
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Table  £ 
Nomenclature 


concentration  Jmoles/ccj 

2 

diffusion  coefficient  - - 

sec. 

diameter  of  lucite  cylinder  (cm) 
electrons 


Faraday  equivalent  96,  500 

equivalent 

dimensionless  stream  function 

function  of  x  used  in  similarity  transformation 


current  density  - - 

cm 

total  current  (ma) 

number  of  electrons  transferred  in  electrochemical  reaction 
Reynolds  number 
Schmidt  number 

total  pressure -atmospheres 

,  2  calories 

universal  gas  constant - ; — 

gm-mole  -  K 


radius 


o 

absolute  temperature  K 

dummy  integration  variable 

x-component  of  velocity  in  boundary  layer 

x-component  of  velocity  at  outer  edge  of  boundary  layer 

free  stream  velocity 


y-component  of  velocity  in  boundary  layer 

co-ordinate  measured  along  circumferential  direction  from  leading 
edge  of  cylinder 
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y 


co-ordinate  measure  radially  from  surface  of  cylinder 
Sc^^^  ri  -  radial  co-ordinate  scaled  by 


Greek  Symbols 


8 

V 
n 

n 

V 

p 


dimensionless  velocity  gradient  at  wall  - 


5u 

dy 


y=0 


gamma  function 

overvoltage -dimensionless  radial  co-ordinate  scaled  by  Re^^®  -  Re^^® 
angular  velocity  (rpm) 


kinematic  viscosity 
density 

ml 


cm 

sec 


9 


similarity  transformation  variable 
dimensionless  concentration  ratio 


g(x) 

c 


Subscripts  and  Superscripts 

i  individual  species 

m  mixture 

s  conditions  at  electrode  surface 

CO  conditions  in  bulk  electrolyte 
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Table  II 


Calculation  of  Diffusion  Coefficients  for 


K^Fe(CN)^  and  K^Fe{CN)^  in  KCl 


•w  o 

b 

Electrode 

Platinized  Platinum 

Area 

0.  42  cm^ 

Voltage 

0.  600^^  N.  H.  E. 

Electrolyte  - 

1.  0  M  KCl 

O  (rpm) 

i(ma) 

-  ma 

~ — 2 
cm'^ 

2.  1 

1.45 

0.  345 

0.  823 

3.  8 

1.95 

0.450 

1.  074 

6.  1 

2.47 

0.  560 

1.  337 

9.  5 

3.08 

0.680 

1.  623 

12.  5 

3.  54 

0.  790 

1.  885 

16.2 

4.03 

0.900 

2.  148 

20.  3 

4.  51 

1.  020 

2.  434 

^Fe(CNL^  =  5.  19 
0 

in-6  2, 

X  10  cm.  / sec  . 

2.  32 

1.  52 

0.  350 

0.  835 

0,  979 

0.  989 

0.  220 

0.  525 

5.  43 

2.  33 

0.  520 

1.  241 

7.  06 

2.66 

0.  590 

1.  408 

3.48 

1.  87 

0.  425 

1.  014 

2.  30 

1.  52 

0.  350 

0.  835 

5.  54 

2.  35 

0.  520 

1.  241 

11,  10 

3.  33 

0.  730 

1.  742 

14.  10 

3.  7b 

0.  810 

1.  933 

18.  20 

4.  27 

0.  920 

2.  196 

^e(CN)^  =5.11 

X  10  cm.  /sec. 
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Table  II  (coat  ) 


Electrode 

-  Platinized  Platinum 

Area 

-  0.42cm^ 

Voltage 

-  -.  175'^N.  H,  E. 

Electrolyte 

-  1.  0  M.  KCl 

fl(rpm) 

i{ma) 

ma 

cm2 

10.  2 

3.  19 

0.  73 

1.  742 

11.62 

3.41 

0.  77 

1. 838 

14.  90 

3.  86 

0.  87 

2.  076 

17.  75 

4.  21 

0.  94 

2.  243 

24.  30 

4.  93 

.1. 

2.  625 

22.  90 

4.  79 

1.  C.50 

2.  506 

21.  25 

4.  61 

1.  025 

2.  446 

19.  70 

4.  44 

0.  980 

2.  339 

26.  10 

5.11 

1.  130 

2.  697 

28.  60 

5.  35 

1.  180 

2.  816 

27.  50 

5.  24 

1.  150 

2.  745 

9.61 

3.  10 

0.  7i0 

1.  695 

8.  64 

2.  94 

0.  670 

1.  599 

7,  SO 

2.  79 

0.  630 

1.  504 

^Fe{CN)^^  " 

6.  20  X  10  ^cm^/sec 
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Table  HI 


Diffusion  Co.  ficients 


Component 

_D 

2 

Supporting  Electrolyte 

Reference 

K^FeCCN)^ 

5.  0  X  io"^^5L. 

sec 

1.  0  M  KCl 

5 

K^Fc(CN)^ 

A  2 

7  A  1  A*®  cin 

7.  0  X  10  - 

sec 

1.  0  M  KCl 

5 

K^Fe{CN)^ 

4.9x 

sec 

0.  5  M  K^SO, 

2  4 

5 

K^FeCCN)^ 

7.  2  X 

sec 

0.  5  M  K^SO, 

2  4 

5 

K^FeCCN)^ 

5.  8  X 

sec 

2  N  NaOH 

2 

K^Fe(CN)^ 

K^Fe(CN)^ 


From  the  Translating  Electrode 


5.  15  X  10 


,-6  cm 


sec 


6.  20  X  10 


-6  cm 


sec 


1.  0  M  KCl 
1.  0  M  KCl 
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Table  IV 


Amplitude  of  Curren 

t  Fluctuations 

Cylinder  Size 

Speed  (rpm) 

Current  (ma) 

4^  X  100 
1 

1” 

82,  8 

0.  73 

19.  1 

35.  7 

0.45 

16.  7 

12.  9 

0.  28 

14.  3 

52.  2 

0.  52 

19.  2 

17.  1 

0.  32 

18.  8 

20.  4 

0.  35 

14.  3 

3/4" 

10.  9 

0.  89 

11.  2 

61.  5 

0.  66 

13.  6 

22.  2 

0.40 

11.  2 

85.  7 

0.  79 

15.  2 

1/4" 

129.0 

0.98 

?.  1 

50.  8 

0.  61 

3.  3 

22.  2 

0.  40 

5.  5 

35.  3 

0.  51 

5.  3 

1/8" 

131.0 

2.  0 

2.  5 

118.0 

1,  83 

2.  7 

24.  2 

0.  80 

3.  8 

5.  35 

0.  42 

5.  2 

v/here  Ai  =  i  -  i  . 

max  mm 
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B.  Heat  and  Mass  Transfer  From  the  Surface  of  a  Cylinder  with 

Discontinuous  Boundary  Conditions  to  an  Incompressible  Laminar 
Flow  (See  Supplement  No.  1  accompanying  this  report.) 
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Figure  3  Frequency  of  Current  Fluctuations  vs.  Rotational  Speed 
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Figure  5  Current  Density  vs.  Square  Root  of  Rotational  Speed 
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Figure  6  Current  Density  vs.  Square  Root  of  Rotational  Speed 


Figure  7  Current  Density  vs.  Square  Root  of  Rotational  Speed 
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Figure  8  Current  Density  vs.  Square  Root  of  Rotational  Speed 


Figure  9  Current  vs.  Square  Root  of  Rotational  Speed 
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Figure  10  Current  vs.  Square  Root  of  Rotational  Speed 


I  (mo/cm*) 


01 _ _ I 

0  2  4  6  8  10  12 

^1/2 


Figure  11  Current  Density  vs.  Square  Root  of  Rotational  Speed 


Figure  12  Current  Density  vs.  Square  Root  of  Rotational  Speed 
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Figure  13  Current  vs.  Square  Root  of  Rotational  Speed 
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Figure  14  Current  vs.  Square  Root  of  Rotational  Speed 


Figure  15  Current  vs.  Square  Root  of  Rotational  Speed 
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5'igure  16  Current  vs.  Square  Root  of  Rotational  Speed 


